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This  research  deals  with  the  implementation  of  a  linear  filter- 
ing system  capable  of  yielding  the  strains  in  a  solid  body  from  the 
information  on  the  displacement  field  contained  in  the  light  intensity 
trace  of  moire  fringes. 

An  analysis  of  the  photo-optical  device  used  to  produce  and 
record  the  fringes  illustrates  quantitatively  the  filtering  abilities 
of  the  lens-film  combination.  Solutions  to  the  non-linear  sensitoraetric 
characteristics  of  the  film  are  applied  to  moire  data. 

Several  experiments  of  optical  filter  synthesis  that  increase 
substantially  the  quality  of  recorded  patterns  are  reported. 

The  retrieval  of  the  meaningful  phase  modulation  of  the  light 
intensity  trace  from  the  spurious  amplitude  and  phase  variations  due  to 
noise  is  successfully  accomplished  by  the  implementation  of  n\imerical 
filters  in  quadrature.  This  is  used  to  determine  new  bounds  on  auto- 
matic fractional  moire  fringe  interpolation. 

xvi 


The  numerical  differentiation  of  the  empirically  determined 
displacement  curve  is  critically  reviewed  and  a  significant  increase 
in  the  accuracy  of  this  process  is  obtained  by  the  use  of  several 
numerical  methods. 

Necessary  techniques  of  model  girid  engravure  able  to  withstand 
high  temperatures  are  developed. 

The  above  improvements  on  the  moire  method  are  used  to  measure 
thermal  strains  in  stainless  steel  rings,  under  steady  state  thermal 
loading,  up  to  a  maximum  temperature  of  1600  F. 
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CHAPTER  I 
INTRODUCTION 

The  moire  method  has  the  inherent  capability  of  performing 
large  field  strain  determinations  in  such  important  areas  of  solid 
mechanics  as  elasticity,  viscoelasticity,  plasticity  and  fracture, 
under  mechanical,  thermal,  static  and  dynamic  loading  conditions. 

In  some  of  these  fields  the  experimental  stress  analyst  has  a 
choice  of  several  applicable  methods,  in  others,  moire  today  offers 
the  only  possible  solution. 

It  is  unnecessary  to  review  here  the  history  of  the  development 
of  the  moire  method  as  an  experimental  stress  analysis  tool.  Complete 
chronological  listings  of  work  performed  on  this  subject  so'e  available 
[1,  2]. 

Though  in  principle  the  method  can  yield  the  information  which 
we  desire  on  the  strain  field  in  the  body,  not  always  does  this  come 
in  a  convenient  and  direct  form.  A  certain  degree  of  sophistication 
in  the  photo-optical  instrumentation  involved  is  therefore  necessary 
if  the  required  degree  of  sensitivity  and  accuracy  in  the  quantitative 
results  is  to  be  achieved. 

In  this  connection  it  is  important  at  this  early  stage  of  our 
presentation  to  stress  as  fundamental  the  difference  between  the  dis- 
crete and  the  continuous  approach  to  the  information  contained  in  the 
moire  fringes  as  was  explained  by  Sciammarella  [  3] . 


The  foundation  of  i*iat  will  be  said  in  this  entire  work  -will  be 
found  in  the  moire  light  intensity-displacement  law  derived  by 
Sciammarella  [3]  and  tested  by  Ross,  Sciammarella  and  Sturgeon  [h] • 

We  accept  therefore  that  the  displacement  component  in  the 
direction  nornial  to  the  master  grid  at  a  point  in  the  deformed  body 
is  related  in  a  continuous  form  to  the  light  intensity  at  that  point 
by  the  expression 


^^^  "  2tt  ^^   ^°^   1 


(1.1) 


where  the  coordinate  x  is  measured  in  the  direction  normal  to  the 
master  grid. 

It  has  been  shown  [U]  that  by  the  proper  use  of  (1.1)  it  is 
possible  to  increase  the  sensitivity  and  accuracy  of  the  moire  method 
using  coarse  grids  of  300  lines  per  inch  far  beyond  the  bounds  set  by 
the  discrete  theory  [5].  FurthexTuore,  it  was  shown  that  the  process- 
ing of  moire  data  could  be  made  almost  totally  objective  vdien  instru- 
mentation was  designed  to  perform  the  operations  made  by  hand  in  the 
discrete  approach  [6]. 

The  principal  goal  of  this  research  was  to  implement  an  experi- 
mental setup  capable  of  performing  quantitative  strain  determinations 
with  sufficient  accuracy  in  metals  at  high  temperatures  by  the  moire 
method. 

Since  a  careful  systems  analysis  of  the  entire  moire  process 
was  thought  to  be  indispensable  to  achieve  our  purpose,  several  aux- 
iliary tests  under  mechanic  loading  had  to  be  performed  whenever  these 


helped  to  clarify  aspects  of  our  analysis.  Traction  tests  on  bars  and 
compression  tests  on  disks  were  made  for  this  purpose. 

The  optical  synthesis  of  moire  data,  not  previously  attempted, 
seemed  to  offer  interesting  possibilities  to  filter  out  the  inevitable 
photographic  noise  present  in  all  moire  patterns  obtained  under  any- 
thing but  the  ideal  conditions  of  illumination  and  surface  character- 
istics. This  matter  was  therefore  pursued  and  definite  improvements  in 
the  quality  of  the  moire  patterns  were  obtained. 

Previous  to  this  research,  an  electronic  band-pass  filter  able 
to  operate  between  0.C2  cps  and  2000  cps,  had  been  used  to  filter  both 
high  and  low  spurious  spatial  frequencies  from  the  moire  light  intensity 
trace  [6].  By  means  of  an  initial  mismatch  the  frequency  of  the  moire 
signal  before  and  after  deformation  was  made  to  fall  within  the  range  of 
the  instrument. 

However,  the  inclusion  of  the  electronic  filter  caused  large 
oscillations  of  long  duration  on  the  analog  voltage  output  of  the  photo- 
reading  device  every  time  a  negative  was  introduced  or  removed  from  the 
electro-optical  circuit.  The  decay  time  of  these  oscillations,  caused 
by  the  response  of  the  instrument  to  a  step  load,  was  found  to  be  in 
excess  of  one  hour. 

Furthermore,  the  displacement  information  due  to  the  deforma- 
tion of  the  solid  body  had  to  be  obtained  as  the  difference  between 
the  final  and  initial  pattern.  When  strong  mismatch  was  present,  in 
conjunction  with  small  deformations,  this  difference  between  two  large 
quantities  of  approximately  the  same  magnitude  impaired  the  accuracy 
of  the  procedure. 


Since  the  amplitude  variations  were  not  corrected  by  the  elec- 
tronic filter,  a  numerical  normalization  technique  was  necessary  in 
the  digital  program  that  handled  the  moire  data  [6  ] .   The  errors 
involved  in  this  step  were  found  to  be  considerable  and  are  discussed 
fully  in  [6]. 

As  an  alternative  to  the  above  procedure,  numerical  filtering 
with  the  aid  of  a  digital  computer  was  introduced  here  as  a  means  of 
retrieving  the  meaningful  spatial  phase  modulation  of  the  light  inten- 
sity trace  from  the  spurious  variations  of  phase  and  amplitude.  This 
was  possible  by  the  implementation  of  in-<}uadrature  filters  that  not 
only  removed  all  the  unwanted  frequency  components  of  the  signal  with 
no  transient  perturbations,  but  also  produced  an  output  completely 
independent  of  the  amplitude  of  the  moire  light  intensity  trace.  This 
was  instrumental  in  establishing  the  practical  bounds  on  the  fractional 
fringe  interpolation  that  is  possible  by  application  of  the  continuous 
law  (1.1). 

The  differentiation  process  that  enables  us  to  obtain  the 
strains  from  the  displacement  field  was  also  given  a  good  deal  of  atten- 
tion. A  critical  evaluation  of  several  differentiation  procedures  of 
experimental  data  was  made.  This  permitted  a  considerable  increase  in 
the  accuracy  of  the  strain  determination  which  is  perhaps  the  single 
most  undesirable  feature  of  the  moire  method  as  applied  to  experimental 
stress  analysis. 

We  mentioned  at  the  beginning  of  this  introduction  several  pos- 
sible applications  of  the  moire  method  in  solid  mechanics;  none  is  more 
challenging  than  the  measurement  of  strains  in  metals  subjected  to 


thermal  loading.  In  this  case  the  absolute  value  of  the  strains  and 
of  their  difference  across  the  body  are  small,  so  both  the  sensitivity 
and  accuracy  of  the  method  of  measurement  are  factors  of  primary 
importance. 

Therefore,  the  ability  of  moire  to  produce  good  results  under 
these  conditions  is  a  suitable  test  of  excellence  for  the  entire  system 
performance.  Such  results  were  obtained  in  this  research. 

The  extension  of  the  moire  method  to  high  temperatures,  at  vriiich 
conventional  photoengravure  techniques  of  model  line  printing  fail 
because  the  photosensitive  materials  used  burn  off,  required  the  estab- 
lishment of  new  methods  of  grid  manufacture.  Several  such  methods  were 
devised  for  aluminum,  copper,  mild  and  stainless  steel. 

By  the  use  of  these  grids,  moire  fringe  patterns  were  photo- 
graphed up  to  a  maximum  temperature  of  1600 °F. 


CHAPTER  II 
THE  PHOTO-OPTICAL  SYSTEM 

2.1  Introduction 

The  photo-optical  system  used  to  produce  and  record  moire 
patterns,  has  the  ability  to  perform  as  a  linear  low-pass  filter  of 
spatial  frequencies. 

The  system  is  composed  of  the  lens  and  the  film,  both  of  vrfiich 
have  individually  this  ability.  Thus,  they  operate  together  as  two 
linear  filters  in  series. 

The  filtering  effect  can  be  thought  of  as  the  process  by  which 
the  higher  spatial  harmonics  of  the  Fourier  decomposition  of  the  object 
are  unable  to  pass  onto  the  recorded  image.  This  has  a  degrading  effect 
on  the  image,  since  the  infinite  expansion  that  represents  the  object 
in  all  its  detail,  is  interrupted  after  a  given  finite  number  of  terms. 

Moreover,  we  will  see  that  the  effect  of  the  system  is  to  damp, 
with  progressive  vigor  as  frequencies  increase,  the  corresponding  ampli- 
tude coefficients  of  the  Fourier  expansion.  This  produces  an  increasing 
loss  of  contrast  as  the  spatial  wavelengths  get  shorter. 

In  -fcrt^at  follows,  we  will  understand  by  system  analysis  the  proc- 
esses by  which  we  record  passively  the  action  of  a  given  system  on  the 
moire  information. 

We  will  call  system  synthesis  the  processes  by  which  we  act  on 
the  system  to  change  its  frequency  response  characteristics  to  suit 
our  purposes. 


2.2  Optical  Analysis  of  the  System 

Due  to  the  filtering  effect  mentioned  above  lens  and  film 
manufacturers  have  sought  to  devise  means  of  evaluating  the  image  form- 
ing capabilities  of  lenses  and  films.  One  such  technique,  based  on  the 
Modulation  Transfer  Function  (MTF),  offers  special  advantages  for  the 
study  of  the  photo-optical  system  used  to  record  moire  data. 

The  method  consists  of  obtaining  the  image  formed  by  the  photo- 
optical  system  of  an  object  made  up  of  parallel  lines  with  a  sinusoidal 
intensity  distribution.  A  typical  test  object  is  manufactured  by  the 
Eastman  Kodak  Company  and  can  be  purchased  commercially. 

The  MTF  can  be  defined  numerically  as  the  ratio  of  the  contrast 
in  the  image  to  the  contrast  in  the  object  for  the  same  spatial 
frequency. 

We  define  as  is  usual  the  contrast  C  of  the  object  test  target 
under  coherent  or  incoherent  illumination,  as  the  ratio  of  the  ampli- 
tude of  the  sinusoidal  intensity  distribution  to  the  mean  background 
luminance,  then 

P   I  max  -  I  min  ron^ 

I  max  +  I  min  ^'^'"^^ 

vfcere  the  mean  value  is  sufficiently  large  so  the  luminance  is  always 
positive. 

If  we  now  define  the  spatial  frequency  of  the  sinusoidal  test 
object  as  uJ  =  2Tt/X^   where  X  is  the  wavelength,  then  the  MTF  is 
expressed  as 

M(ou)  =  tlnl  (2.2) 


vrfiere  C  (ua)  and  C  (uu)  are  the  contrast  of  the  image  and  the  object  at 
frequency  tu. 

The  MTF  is  then  a  measure  of  the  relative  degradation  in  ampli- 
tude of  the  image  of  a  sinusoidal  object. 

Under  incoherent  illumination,  a  condition  common  in  moire  xrork, 
the  optical  response  of  the  lens  and  the  film  may  be  regarded  as  two 
linear  operations  in  series.  This  requires,  however,  that  the  exposure 
level  of  the  photograph  be  adjusted  so  that  the  entire  amplitude  of  the 
sinusoidal  input  will  lie  within  the  straight  line  portion  of  the  H  &  D 
characteristic  curve  of  the  film. 

Under  the  above  conditions  of  linearity,  the  MTF  of  the  photo- 
optical  system  can  be  expressed  as 


M 


:^(a))  =  Mj^(u;)  .  Mj,(a))  (2.3) 

where  M  (uu),  M^  (cd)  and  M„(uj)  are  the  transfer  functions  of  the  photo- 
optical  system,  the  lens  and  the  film,  in  that  order. 

The  MTF  technique  offers  the  following  advantages  for  moire 
work: 

a)  The  test  targets  used  in  the  method  resemble  exactly  the 
actual  objects,  in  fact,  the  sinusoidal  test  targets  can  be  considered 
as  perfect,  noiseless,  undlstorted  moire  fringes. 

b)  The  performance  of  the  system  is  evaluated  over  a  large  area 
in  the  object  plane. 

c)  The  MTF  can  not  only  be  used  to  predict  the  response  of  the 
entire  system,  but  can  also  be  used  to  examine  each  element  in  the 


;om- 


photo-optical  process  to  determine  what  degradation  of  the  image  is 
introduced  at  each  stage. 

d)  The  MTF  of  a  system  can  be  obtained  from  those  of  its  cc 
ponent  elements  by  simply  multiplying  the  corresponding  MTF  curves, 
ordinate  by  ordinate,  for  each  abscissa  value. 

e)  The  knowledge  of  the  MTF  curves  of  the  system,  as  will  be 
explained  later,  enables  us  to  determine  the  largest  maximum  displace- 
ment that  can  be  measured  by  the  system. 

f)  For  the  critical  work  of  reproduction  of  moire  line  masters, 
the  MTF  offers  again  quantitative  information  regarding  the  contrast 
that  can  be  expected  between  dark  and  clear  lines  in  the  copy. 

g)  Finally,  the  MTF  curve  adequately  represents  the  optical- 
filtering  capabilities  of  the  system. 

Notwithstanding  the  above,  Hempenius  [7]  has  found  that  photo- 
graphic systems  analysis  based  only  on  the  MTF  characteristics  may  be 
in  some  cases  significantly  in  error  and  that  the  effect  of  grain  should 
be  taken  into  account.  This  is  not  believed  to  be  the  case  in  moire 
fringe  photography. 

2.3  Modulation  Transfer  Function  of  Lenses 

Camera  lenses  degrade  the  arriving  light  signal  by  aberration 
and  by  diffraction. 

We  will  first  consider  the  diffraction  effects,  assuming  for 
this  purpose  that  we  possess  an  ideal  aberration-free  lens. 

Between  the  spatial  wavelengths  large  enough  to  be  passed  with- 
out attenuation  of  amplitude  and  the  one  that  represents  the  limit  of 
resolution,  the  MTF  of  an  ideal  lens  varies  as  shown  in  Figure  la. 
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We  vdll  now  turn  to  the  aberration  losses,  which  depend  on  the 
lens  construction.  Photographic  lenses  must  be  made  to  operate  at 
various  aperture  settings  and  the  best  any  manufacturer  can  do  is  offer 
several  different  compromises  between  the  various  types  and  amounts  of 
aberrations:  distortion,  astigmatism,  curvature  of  the  field,  etc. 
However,  the  optimum  chromatic  correction  available  should  be  preferred 
for  both  moire  master  and  fringe  photography. 

Regarding  the  other  aberration  losses,  unfortunately  no  recom- 
mendations of  such  a  general  nature  can  be  given.  Kelly  [8]  reports 
tests  on  a  50  mm  f/2.8  photographic  lens  that  show  that  geometrical 
aberrations  are  responsible  for  most  degradation  at  f/2,8  v±iile  at  f/8 
most  degradation  is  due  to  diffraction.  Hence  the  aperture  setting  of 
the  iris  has  an  important  effect  on  the  MTF. 

This  situation  obliges  us  to  resort  to  actual  testing  of  the 
lens  we  intend  to  use  for  a  future  optical  filtering  application. 
Modification  to  the  curve  of  Figure  la  due  to  aberration  must  be  deter- 
mined experimentally  for  the  magnification  and  aperture  setting  we 
intend  to  use  if  this  information  is  not  supplied  by  the  lens  manufac- 
turer . 

A  test  of  this  nature  is  described  in  Section  2.5  in  detail. 

2.I4.  Modulation  Transfer  Function  of  Films 

The  Modulation  Transfer  Function  of  a  film  decreases  with 
increasing  spatial  frequency  and  resembles,  in  a  general  way,  that  of 
a  lens.  The  damping  of  high  frequency  amplitude  coefficients  reflects 
the  effect  of  the  diffusion  of  light  within  the  emulsion  on  the  micro- 
structure  of  the  image. 
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However,  all  films  do  not  possess  the  same  limiting  curve  of 
maximum  resolution  as  is  the  case  with  lenses.  Films  with  a  MTF 
approximately  equal  to  one  over  a  broad  range  of  spatial  frequencies 
are  manufactured  for  spectroscopy  and  document  reproduction. 

The  MTF  for  the  Kodak  Royal  Ortho  film  used  in  the  photo- 
optical  test  described  in  Section  2.5  can  be  seen  in  Figure  lb. 

Complete  information  on  the  MTF  of  numerous  photographic  films 
and  plates  is  available  from  the  manufacturers  [9]. 

2.5  Experimental  Determination  of  the  Modulation 
Transfer  Function  of  the  Photo-Optical  Syst'em 

Tests  were  conducted  to  detennine  the  MTF  of  a  12  in.,  f/10, 
Alphax  copying  lens,  and  of  this  lens  together  with  Kodak  Royal  Ortho 
film. 

The  lens  and  film  were  used  later  in  the  thermal  test  of 
Chapter  IV.  The  magnification  tested  was  one  to  one  and  the  aperture 
setting  was  f/10. 

Since  the  MTF  will  vary  across  the  image  plane  of  a  photo- 
graphic object,  from  its  highest  value  near  the  optical  axis  to  a  lower 
value  near  the  edges,  test  targets  were  placed  in  the  center  and  at  one 
edge  of  the  object  plane. 

The  test  targets  manufactured  by  the  Eastman  Kodak  Company  vary 
sinusoidally  in  transmittance,  with  harmonic  distortion  of  less  than 
2^,  and  range  in  spatial  frequency  from  3/8  to  ii2  cycles  per  millimeter. 
The  modulation  of  transmittance  is  approximately  6S%   for  the  target 
placed  in  the  center  of  the  field  and  30^  for  the  one  placed  at  the  edge. 
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The  targets  also  possess  a  grey  scale  for  calibration  purposes.  The 
optimum  performance  of  the  lens  was  found  to  be  restricted  to  a  circle 
with  center  at  the  optical  axis  and  radius  0.16  the  focal  distance. 

Since  the  focussing  of  the  camera  has  a  great  influence  on  the 
MTF  the  exact  focus  will  provide  the  highest  limit  of  resolution.  For 
every  position  other  than  the  true  focus  the  lens  performs  with  a  differ- 
ent  MTF.  The  further  the  inage  plane  is  removed  from  the  focal  plane, 
the  poorer  the  resolution  becomes  [10]. 

This  immediately  suggests  a  way  of  filtering  out  excessive  low 
frequency  passed  by  the  lens  by  slightly  defocussing  the  optical  system. 

This  also  obliges  us  to  perform  the  lens  test  for  the  same  posi- 
tion of  the  image  plane  with  respect  to  the  true  focal  plane  as  that 
at  which  we  will  operate  during  the  actual  running  of  the  moire  test. 
This  is  not  an  easy  matter  since  both  tests  may  be  performed  at  differ- 
ent locations  and  with  different  optical  equipment  except  for  the  lens 
tested. 

The  recommended  technique  [11]  of  taking  several  photographs  at 
positions  close  to  the  true  focus  and  selecting  the  optimum  between  them 
is  too  lengthy  for  most  experimental  stress  analysis  applications  where 
the  object  may  have  to  be  moved  several  times  during  the  test. 

We  have  used,  to  select  the  optimum  focussing  in  both  the  lens 
and  thermal  tests,  the  visual  inspection  with  the  aid  of  a  60X  micro- 
scope of  the  image  of  12  lines  per  millimeter  produced  on  a  finely 
ground  glass  manufactured  by  us  for  this  purpose. 

To  record  the  image  we  selected  Kodak  Royal  Ortho  film  with  an 
ASA  speed  of  i|00  because  of  its  known  MTF,  H  &  D  characteristic,  and 
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proven  ability  to  record  moire  patterns  at  elevated  temperatures. 
Care  was  exercised  to  expose  the  sinusoidal  test  patterns  entirely  in 
the  straight  line  portion  of  the  sensitometric  curve. 

The  MTF  of  the  film  is  represented  in  Figure  Ic.   The  MTF  of 
the  lens-film  combination  was  obtained  by  scanning  both  the  test  object 
and  the  negative  with  a  Joyce  and  Loebl  microdensitometer .  By  scanning 
the  grey  scales  of  the  target,  a  calibration  curve  can  be  obtained 
relating  the  densitometer  reading  to  the  corresponding  density  [6]. 

The  scans  of  the  sinusoidal  patterns  were  then  analysed  to 
determine  the  average  values  of  the  maximum  and  minimum  densities  for 
a  given  frequency. 

The  modulation  was  then  obtained  by  the  use  of  the  formulas  [12] 

AD  =  log  T  min  (uu)  -  log  T  max  (u))  (2.U) 

C(aO  =  T  max  (u))  -  T  min  (u;)  ,      ^. 

^    '        T  max  (ud)  +  T  min  (uj)  ^'^'^> 

The  plot  of  the  modulation  transfer  function  M  (u))  correspond- 
ing to  the  lens-film  combination  that  was  obtained  from  the  experiment 
can  be  seen  in  Figure  Id. 

According  to  (2.3),  if  we  have  measured  M  (aO  and  we  are  given 

o 

MpC'u),  the  MTF  of  the  lens  alone  can  be  obtained.  The  results  have 
been  plotted  in  Figure  le. 

The  product  of  the  MTF  of  the  ideal  lens  and  the  film  differs 
considerably  from  the  measured  values  of  the  combination.  This  is  par- 
ticularly so  as  the  wavelengths  get  shorter.  The  vigor  of  this  damping 


Ill 


is  undoubtedly  due  to  incorrect  focussing  and  aberrations  rather  than 
to  the  latter  alone.  This  emphasizes  the  strong  dependence  of  the 
transfer  characteristics  of  the  lens  on  the  correct  focussing  [10]. 

Figure  If  shows  the  results  of  testing  the  lens-film  combina- 
tion at  one  edge  of  the  image  plane.  Comparison  of  Figure  Id  with 
Figure  If  shows  how  the  system  performs  better  for  smaller  angular 
apertures.  The  damping  of  the  very  long  wavelengths  indicates  loss  of 
focus  by  lack  of  parallelism  of  the  image  plane  with  respect  to  the 
object  plane. 

2.6  Light  Intensity  Trace  of  Moire  Fringes 

The  light  intensity  variation  along  a  line  in  a  moire  pattern 
given  by 

I(x)  =  Iq  +  !]_  cos  ^n  p(x)  (2.6) 

can  be  obtained  from  the  photographic  negative  used  to  record  the 
pattern.  A  photo-reading  device  as  the  one  presented  in  [h]   can  be 
used  to  obtain  the  trace  of  the  transmittance  of  the  negative,  or 
a  scanning  microdensitometer  can  be  employed  to  record  the  density  varia- 
tions of  the  film. 

The  incident  exposure  on  the  film  can  then  be  retrieved  from  the 
transmittance  or  from  the  density  trace  by  performing  the  operations 

_1_ 

E  (x)  =  I(x).t  =  T„(x)^^"  (2.7) 

n  n 


1 


E  (x)  =  I(x)-t  = 
n  ' 


Log'^  D^(x) 


^  (2.8) 


as  the  case  may  be, 
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It  has  been  the  custom  in  moire  work  to  regard  the  negative, 
obtained  by  photographing  the  fringes,  as  a  suitable  representation  of 
the  original  exposure  pattern.  The  clear  fringes  in  the  model  appeared 
as  dark  fringes  in  the  negative  and  vice  versa,  but  no  consequence  was 
attached  to  this  since  the  curve  of  displacements  was  plotted  with 
respect  to  an  arbitrary  origin. 

A  look  at  the  characteristic  curve  of  a  film  shows  immediately 
that  the  output  of  a  film  to  incident  exposure  is  transmittance  or 
density  of  silver  deposit.  Furthermore,  this  relationship  is  non- 
linear and  independent  of  the  spatial  frequency  of  the  input.  There- 
fore, vdiat  is  recorded  on  film  of  an  arriving  spatial  sinusoidal  varia- 
tion of  intensity  is  a  harmonically  distorted  waveform  of  spatial  den- 
sity or  transmittance,  Figure  2. 

The  amount  of  this  distortion  is  illustrated  in  Figure  3  vdiere 
exposure  rather  than  the  logarithm  of  the  exposure  has  been  plotted  in 
the  abscissas.  In  Figure  h   we  show  the  strong  distortion  of  an  actual 
transmittance  trace  of  a  sinusoidal  light  intensity  variation  corre- 
sponding to  a  moire  pattern  of  fringes. 

Since  the  application  of  the  continuous  light  intensity  dis- 
placement law  makes  use  of  the  light  intensity  of  all  portions  of  the 
spatial  waveform,  the  undistorted  reproduction  of  the  signal  is  manda- 
tory and  Equations  (2.7)  or  (2.8)  should  be  applied  always. 

It  will  be  seen  in  Section  2.7  that  it  is  possible  to  perform 
photographically  these  operations. 

It  is  interesting  to  note  that  there  is  no  space  shift  of  the 
maximums  and  minlraums  of  the  intensity  curve  with  respect  to  those  of 
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the  transmittance  curve.  Therefore,  the  usual  method  of  hand  analysis, 
that  only  takes  these  points  into  account  for  the  plot  of  the  displace- 
ment curve,  will  be  free  of  the  errors  due  to  photographic  reproduction 
discussed  here. 

The  continuous  trace  of  incident  exposure  will  contain,  in  addi- 
tion to  the  phase  modulated  spatial  curve  containing  the  information  on 
the  displacement  in  the  body  given  by  (2.6),  spurious  amplitude  and  phase 
variations  due  to  noise.   These  disturbances  are  visible  in  Figure  h. 

¥e  wish  to  extract  (2.6j  from  the  noise  in  which  it  is  submerged 
before  proceeding  with  the  determination  of  displacements  to  increase 
the  accuracy  of  the  procedure. 

The  amplitude  variations  are  caused  by  the  non-uniform  illumina- 
tion of  the  object,  by  changes  in  the  reflectance  of  the  model  surface, 
and  by  the  progressive  vigor  with  which  the  amplitude  coefficients  of 
the  Fourier  decomposition  of  the  image  are  damped  by  the  filtering 
action  of  the  photo-optical  system. 

The  more  frequent  causes  of  short  wavelength  variations  are  the 
image  of  the  master  and  model  grid  lines  on  the  moire  fringe,  the 
scratches  and  pits  in  the  model  surface,  and  the  dust  that  can  settle 
over  the  model  during  testing.  The  long  wavelength  variations,  on  the 
other  hand,  are  often  due  to  non-uniform  illumination  and  surface 
reflectance  or  transmission. 

Since  these  perturbations  in  the  trace  are  of  both  longer  and 
shorter  wavelengths  than  those  of  the  useful  signal,  a  band-pass  filter 
is  necessary  to  eliminate  them. 
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Optical  band-pass  filters  cannot  be  easily  synthesized;  however, 
the  photo-optical  system  used  to  produce  and  record  the  moire  fringes 
has  inherent  low-pass  filter  capabilities  of  which  we  can  take  advantage. 
This  will  allow  us  to  implement,  optically  at  least,  one  part  of  our 
filter  requirements. 

Both  long  wavelength  and  amplitude  perturbations  may  be  removed 
by  other  than  optical  means,  and,  to  do  this,  digital  In-quadrature 
filters  have  been  used  successfully  as  will  be  seen  in  Chapter  III. 

It  remains  for  us  to  develop  in  this  chapter  the  techniques  to 
photographically  correct  the  wave  distortion  due  to  the  non-linear  film 
response  and  to  substantially  reduce  optically  the  noise  content  of  the 
moire  fringe  light  intensity  trace. 

2.7  Correction  of  Spatial  Waveform  Distortion 

The  so-called  straight  line  portion  of  the  characteristic  curve 
of  the  film  is  a  linear  relationship  between  the  density  and  the  logar- 
ithm of  the  exposure.  The  logarithm  of  a  sinusoidal  exposure  variation 
is  shown  in  Figure  5. 

Therefore,  a  sinusoidal  exposure  within  the  straight  line  portion 
of  the  sensitometric  curve  will  produce  a  spatial  density  variation  lin- 
early proportional  to  the  logarithm  of  the  arriving  sinusoidal  spatial 
waveform.  However,  if  the  exposure  cannot  be  kept  entirely  within  this 
straight  section  of  the  sensitometric  curve  and  the  toe  of  this  curve 
is  used  to  record  the  arriving  signal,  further  distortions  are  produced 
as  shown  in  Figure  6. 
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If  we  wish  to  eliminate  photographically  the  spatial  harmonic 
distortion  due  to  the  non-linear  response  of  the  negative  film,  the 
following  two  processes  have  been  proven  effective. 

a)  Negative-positive  film  combination 

A  combination  of  negative  and  positive  films  can  be  chosen  such 
that  the  products  of  their  gamma  coefficients  is  equal  to  one  [13]. 

For  the  straight  line  portion  of  the  characteristic  curve  of 
the  negative,  the  following  holds 

D  =  V  log  E  +  C^  (2.9) 

n    n   ^  n    1  \     ■' > 

where  D  is  the  density  of  the  developed  film,  E  is  the  exposure, 
Y  the  slope  of  the  sensitometric  curve  which  in  general  will  be  posi- 
tive and  larger  than  one,  and  C-.  is  a  constant.  The  subscript  n  refers 
to  the  negative. 

The  transmission  of  the  negative  is  given  by 

-V 

T^  =  ^2^Xi^  ^2-^°^ 

where  C_  is  a  constant.  If  we  print  this  negative  on  a  positive  film 

with  a  V  such  that 
P 

y^Vp=l  (2.11) 

we  obtain 

T  =CoE  e"P=CoE  (2.12) 

p    2  n         2  n 

where  T  is  the  transmission  of  the  positive  and  C„  is  a  constant, 
p  ^2 

Thus  a  film  linearly  proportional  in  transmission  to  the  exposure  of 
the  object  has  been  obtained. 
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This  process  cannot  always  be  applied  in  moire  work  because  of 

the  diversity  of  films  that  are  used  as  negatives  depending  on  the 
applications  and  the  very  few  pairs  of  films  -v^iose  product  of  gammas 
is  equal  to  one. 

In  the  thermal  test  reported  in  Chapter  IV,  for  example,  six 
different  films  were  used  to  record  the  fringe  patterns,  since  varia- 
tions in  the  ASA  speed  became  necessary  as  temperatures  increased. 

Since  the  value  of  y   depends  on  the  developing  time,  some  suit- 
able pairs  of  films  can  be  obtained  if  the  stringent  requirements  on 
processing  can  be  enforced. 

b)  Non-linear  density  wedge  method 

In  this  method  we  make  use  of  a  microdensitometer  as  the  Mark 
III  C  manufactured  by  Joyce  Loebl  and  Company  of  Gatehead,  England. 
This  machine  measures  density  on  the  negative  by  balancing  it  with 
a  known  density  located  on  a  linear  density  wedge. 

The  principle  of  operation  is  based  on  a  true  double-beam  light 
system  in  which  two  beams  from  a  single  light  source  are  switched  alter- 
nately to  a  single  photo-multiplier. 

One  beam  of  light  is  passed  through  a  point  on  the  film,  the 
other  through  a  calibrated  reference  x^edge  of  linear  density  variation. 
If  the  two  beams  arriving  at  the  photo-amplifier  are  of  different  inten- 
sity, a  servo  motor  is  activated  that  causes  an  optical  attenuator  to 
reduce  this  difference  to  zero.  In  this  way  a  continuously  null  balanc- 
ing system  is  obtained  in  which  the  position  of  the  optical  alternator 
is  made  to  record  the  density  at  any  point  of  the  specimen. 
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Under  the  above  mentioned  conditions  the  machine  will  perform 
linear  density  measurements.  However,  since  we  are  interested  in 
retrieving  the  exposure  sinusoidal  variation,  the  density  wedge  of  the 
machine  can  be  modified  to  achieve  this. 

Figure  7  shows  the  density  trace  of  a  sinusoidal  exposure,  the 
sensitometric  characteristic  curves  of  the  linear  wedge  and  of  a  non- 
linear wedge  with  the  corresponding  output  functions. 

The  linear  wedge  produces  a  tme   density  measurement  at  each 
point  which  results  in  a  distorted  exposure  curve.  The  non-linear  wedge 
will  produce  an  output  identical  to  the  exposure  input,  if  the  sensito- 
metric curves  of  the  film  and  wedge  are  identical. 

To  accomplish  this  practically,  we  have  printed  on  the  type  of 
film  being  used  as  a  negative,  the  appropriate  linear  density  wedge. 
The  developing  process  -was  identical  for  both  the  negative  and  the  non- 
linear wedge  so  the  resulting  y   coefficient  of  both  sensitometric 
curves  was  the  same.  The  results  obtained  by  this  method  are  shown 
in  Figure  8. 

The  linear  wedges  photographed  on  film  were  placed  between  two 
glass  plates  held  together  with  Canada  Balsm.  The  resulting  non-linear 
wedges  were  placed  in  the  travelling  carriage  of  the_,microdensitometer, 
which  then  performed  linear  measurements  of  exposure  for  this  particular 
film. 

The  advantage  of  this  procedure  over  the  negative-positive  film 
combination  is  that  the  appropriate  non-linear  wedge  characteristic 
curve  can  always  be  obtained  for  any  negative.  Both  share  the  dis- 
advantage of  requiring  strict  film  processing  control. 
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2.8  Filtering  of  a  Moire  Signal 

Several  successful  optical  filtering  operations  have  been 
reported  in  the  communication  sciences  and  optical  literature  [lU,  1$, 
16],  however,  moire  work  has  not  previously  benefitted  from  these  use- 
ful techniques. 

¥e  will  describe  in  -v^at  follows  several  experiments  in  which 
the  quality  of  moire  fringes  has  been  substantially  Increased  by 
optical  filtering. 

As  we  said  in  Section  2.6,  an  optical  filter  of  the  low-pass 
type  can  be  readily  Implemented,  in  fact,  the  tests  reported  in  Sec- 
tion 2.5  indicate  that  a  perfect  lens  has  high  frequency  filtering 
ability.  However,  the  cut-off  frequency  of  a  real  lens-film  combina- 
tion may  not  be  sufficiently  short  to  eliminate  the  undesirable  disturb- 
ances of  the  light  intensity  trace  of  the  moire  fringes.  Therefore, 
additional  modifications  of  the  transfer  characteristics  of  the  optical 
system  are  necessary. 

The  purpose  of  this  section  is  to  establish  the  practical  means 
of  achieving  these  modifications  in  moire  work.  It  will  be  seen  that 
strong  disturbances  are  not  always  completely  eliminated  optically  but 
that  the  noise  level  is  substantially  reduced  in  all  cases. 

We  should  keep  in  mind  during  this  discussion  the  structure  of 
the  image  of  a  moire  fringe  shown  in  Figure  9a.  The  phase  variations  of 
the  average  light  intensity  trace  indicated  in  this  figure  carry  the 
information  regarding  the  displacements  in  the  body.  Our  purpose  is  to 
obtain  this  average  by  optically  erasing  the  image  of  the  individual 
grid  lines  and  the  random  high  frequency  variations  due  to  other  sources 
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of  noise.  Failure  to  do  this  will  result  in  an  intolerable  amount  of 
noise  in  the  analog  voltage  output  of  our  photo-reader  that  is  to  be 
automatically  processed  to  yield  the  displacements  and  strains. 

It  is  necessary  to  distinguish  in  our  treatment  between  the 
coherent  and  incoherent  illumination  over  the  object  plane.  In  prac- 
tice, conditions  of  illumination  are  such  that  some  degree  of  coherence 
is  always  present,  but,  for  simplicity,  only  the  extreme  cases  mentioned 
above  will  be  assumed  here  to  exist. 

a)  Coherent  illumination 

Coherent  light  waves  add  in  amplitude  and  linear  relations  hold 
between  the  amplitudes  at  different  points  in  the  system.  Furthermore, 
these  systems  possess  the  property  that  a  Fourier  transform  relation 
exists  between  the  front  and  rear  focal  plane  of  a  lens. 

Under  these  conditions  the  amplitude  distribution  A  over  the 
output  or  image  plane  is  given  by  Figure  10a. 

eo 

A,(x3,y3)  =  n  \J\>y^)\(-r\>  ^3-^1)^1  ^^1  (2-^^) 

vrfiere  A,   is  the  amplitude  distribution  over  the  input  or  object  plane 
and  h.  is  the  amplitude  response  to  a  point  source.  In  the  spatial 
frequency  domain  the  above  equation  is  equivalent  to 

I  (u>  ,u)  )  =  H.  (u;  ,u)  )  A,  (w  ,cu  )  (2.1ii) 

o  x'  y'    A^  x'  y'      in  x*  y' 

where  A  ,H, ,!..      are  the  Fourier  transforms  of  A  ,  H.  and  A,  , 
o'  A^  m  o'  A     in 

respectively. 
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H.  is  therefore  the  amplitude  transfer  function  of  the  optical 
system  that  operates  on  the  object  amplitude  spectrum  to  give  the  image 
amplitude  spectrum. 

Due  to  the  Fourier  transform  relation  that  holds  between  the 
front  and  rear  focal  plane  of  a  lens.  Figure  10a,  the  function  H.  is 
displayed  at  Fo,  thus  affording  a  means  of  controlling  the  Fourier  com- 
ponents that  make  up  the  image  by  inserting  an  appropriate  mask  in 
this  plane. 

The  mask  can  in  general  be  a  complex  transmittance  M  =  |M|e  , 
thus  affording  control  over  the  amplitude  by  varying  the  optical  den- 
sity of  the  mask  and  over  the  phase  by  varying  the  thickness.  One  such 
complex  filter  function  has  been  implemented  recently  by  Vander  Lugt  [17] 

In  dynamic  and  high  temperature  moire  investigations  it  is  often 
necessary  to  produce  a  fringe  pattern  between  the  image  of  the  model's 
grid  and  a  master  grid  placed  in  the  back  of  a  photographic  camera. 
The  moire  pattern  thus  recorded  on  film  contains  not  only  the  image 
of  the  fringes  but  also  the  contact  copy  of  the  line  master.  These 
lines  produce  an  undesired  high  frequency  disturbance  in  the  light 
intensity  trace  of  the  fringes,  Figure  9b. 

When  the  master  is  in  the  object  plane,  a  well-focused  photo- 
graphic camera  will  also  record  the  image  of  the  individual  lines 
together  with  the  moire  fringes,  if  the  aperture  is  not  sufficiently 
reduced  to  perform  a  type  of  filtering  that  will  be  discussed  in  (b) . 
In  this  case  again,  the  images  of  the  grating  lines  are  a  source  of 
noise,  Figure  9a. 
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To  remove  this  disturbance,  the  negative  of  the  moire  pattern 
was  inserted  in  plane  P  of  Figure  10a.  The  spectrum  of  the  image  of 
the  grating  lines  was  displayed  at  Pp  and  a  completely  opaque  mask 
wide  enough  to  interrupt  the  central  order  of  the  Fourier  decomposition 
was  placed  at  Pp.  The  remaining  orders  were  collected  by  the  lens  L^ 
to  produce  an  image  almost  completely  free  of  the  grating  lines  at  P^. 

A  50X  magnification  of  the  appearance  of  the  micro strricture  of 
the  image  befdre  the  bias  removal  is  shown  in  the  upper  insert  of  Fig- 
ure 11,  v^ile  the  lower  insert  corresponds  to  a  50X  magnification  after 
filtering. 

The  strong  damping  is  due  to  the  fact  that  the  film,  as  the  eye, 
is  unable  to  register  negative  intensities  and  records  only  a  pattern 
proportional  to  the  square  of  the  amplitude  of  the  arriving  disturbance. 

b)  Incoherent  illumination 

Incoherent  light  waves  add  in  intensity  and  linear  relations 
hold  between  the  intensities  at  different  points  in  the  system. 

The  intensity  distribution  over  the  image  or  output  plane 
I  (X2>y2)  is  given  by  Figure  10c. 

00 

I  (xg^yg)  =  a  Iii^»yi)N:(^2-'^'  y^-y^)dx^dy^  (2.15) 

where  I.  (x.,y^)  is  the  intensity  distribution  over  the  object  or 
input  plane  and  h^  is  the  intensity  response  to  a  point  source.  In  the 
spatial  frequency  domain  the  above  equation  is  equivalent  to 

T  (cu  ,u)  )  =  T,  (o)  ,<J0  )!L(u)  ,u)  )  (2.16) 

o  x'  y^    in  x'  y^  i  x'  y' 
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where  T  ,  T.  ,  iL  are  the  Fourier  transforms  of  I  ,  I.   and  H,., 
0*  in'  \[  o'  m     I' 

respectively. 

H-.  is  then  the  intensity  transfer  function  that  operates  on  the 
object  intensity  spectrum  to  yield  the  image  intensity  spectmm. 

It  can  be  shown  [18]  that 


/      UJ         (1) 


"l^VV  =  it)  ^^^   ^^^'^^  ^  "^  -  T'  e  -  ■/  d«  dP       (2.17) 

where  f(Q',  3)  is  the  complex  amplitude  distribution  over  the  surface 
of  the  aperture  and  f  (o", 3)  its  complex  conjugate.  The  angular 
coordinates  o/   and  P  are  shown  in  Figure  10b.  C  is  a  constant  and 

0 

K  =  2n/X^  where  X  is  the  wavelength  of  the  radiation. 

Therefore,  control  over  the  amplitude  and  phase  of  f(a,3) 
affords  a  means  of  achieving  various  filtering  effects. 

For  a  perfect  lens  f(Q',3)  is  real  and  equal  to  one  over  the 
entire  aperture  and  zero  elsewhere.  The  application  of  (2.10)  to  a 
perfect  lens  results  in  the  curve  of  Figure  la  and  indicates  that  a 
perfect  lens  is  limited  in  resolution  by  its  aperture.  The  limit  of 
resolution  being  equal  to  d/X .  This  immediately  suggests  a  simple 
way  of  damping  the  high  frequencies  present  in  the  object  by  reducing 
the  size  of  the  iris  of  the  optical  instrument. 

Application  of  this  filter  effect  to  moire  fringes  produced  the 
results  shown  in  Figure  12  vrtiere  the  microdensitometer  traces  of  two 
fringes  corresponding  to  the  same  traction  test  are  presented.  The 
master  grid  was  in  the  object  plane  when  the  photograph  was  taken. 
The  upper  fringe  corresponds  to  an  aperture  setting  of  f/10  and  the 
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lower  to  an  aperture  of  f/6ii.  The  damping  of  the  300  line  per  inch 
of  the  master  is  clearly  visible. 

This  method  is  not  always  practical  in  moire  work  because 
exposure  times  become  very  long  for  the  small  apertures  necessary  to 
filter  out  coarse  grids  at  one  to  one  magnification.  However,  since 
smaller  apertures  will  always  filter  some  high  frequency  noise,  though 
not  perhaps  the  grid  lines,  they  should  be  used  -vriienever  possible. 

We  next  will  take  advantage  of  the  strong  influence  of  defocus- 
sing  on  the  transfer  function  of  a  lens  which  was  already  evident  in 
the  tests  of  Section  2. 5. 

When  we  move  the  image  plane  from  the  focal  plane  we  leave  the 
Frauenhofer  diffraction  region  and  enter  the  Fresnel  zone  vrtnere  the 
diffraction  pattern  at  the  output  plane  will  no  longer  be  the  Fourier 
transform  of  the  aperture.  However,  Cheatham  and  Kohlenberg  [18], 
using  the  approximations  of  geometrical  optics,  obtain  as  the  unit 
intensity  response  of  the  defocussed  system 


/ 


l(vf  ,  ^  |a| 
-  —  for  r  <  - — ' 
a  u  /         -  V 


h(x,y)  =  h(r)  =  \  (2.18) 

0      f  or  r  >  i^ 


where  I  is  the  intensity  of  the  point  source  per  solid  angle,  u  the 
object  distance,  v  the  image  distance  and  a  the  defocal  length. 
Since 

"»        -i(iu  +u)  ) 
H(u)  u,  )  =  If   h(x,y)e    ^  ^^  dx  dy  (2.19) 
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then 


H(a)r)  =  2tt  I'^'   ^ 


aLuur 


ui    aLujr 

V 


(2.20) 


where 


Vx+y   ,     (jur=v<^  + 


2 
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According  to  (2.13)  an  increase  in  the  defocal  distance  -will 
close  the  low-pass  filter  characteristic  of  the  system.  The  practical 
significance  of  this  fact  can  be  seen  in  Figure  13  where  the  micro- 
densitometer  traces  corresponding  to  the  same  fringe  of  a  traction  test 
has  been  photographed  at  different  focal  settings.  The  upper  fringe  of 
the  picture  was  taken  at  the  correct  focus,  the  lower  with  the  image 
plane  removed  a  =  0.0625  in.  from  the  focal  plane  of  the  lens. 

Further  experimentation  was  conducted  to  increase  the  quality  of 
the  image  of  periodic  line  patterns,  whether  they  were  moire  fringes  or 
grid  lines,  under  incoherent  illumination  and  submerged  in  strong  random 
noise. 

In  this  connection  there  exists  an  extensive  literature  in  the 
field  of  optics  dealing  with  the  calculation  of  the  intensity  transfer 
function  H-Cai  ,u)  )  of  a  lens  when  aberrations  are  present.  Steel  [19] 
and  0'Neill[20]  have  calculated  such  functions  for  an  annular  aperture 
and  Francon  [^1]  reports  work  done  by  Sayanagui  on  a  circular  aperture 
covered  by  small  randomly  distributed  transparent  elements. 

The  family  of  modulation  transfer  functions  obtained  by  the 
first  two  authors,  for  different  ratios  of  interior  to  exterior  diameter 
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of  the  annulus,  show  how  the  transfer  characteristics  of  a  lens  can 
be  modified  to  obtain  essentially  different  filtering  effects. 

As  the  ratio  of  the  diameters  approaches  one,  the  transfer 
function  develops  narrow  peaks  in  the  vicinity  of  uu  =  0  and  u)  =  2P  , 
Inhere  3  =  kd/2r,  k  =  2n/X  and  d  is  the  diameter  of  the  aperture. 

Since  this  effect  would  be  useful  to  eliminate  random  noise 
from  a  periodic  line  pattern  which  had  a  spatial  frequency  equal  to  23  , 
experiments  were  conducted  with  annular  apertures  of  the  type  described 
above.  However,  the  loss  of  illumination  for  ratios  of  internal  to 
external  diameter  approaching  unity  was  so  severe  that  photography  of 
the  filtered  pattern  became  impossible. 

Turning  now  to  the  work  done  by  Sayanagui,  this  investigator 
produced  the  random  dots  of  vacuum  deposition  on  a  flat  piece  of  glass 
with  parallel  faces  v*iich  he  then  placed  before  the  lens.   By  this  means 
he  produced  a  filter  that  could  eliminate  in  the  image  plane  P,  of  an 
incoherent  light  system.  Figure  10c,  almost  all  traces  of  the  screen 
mesh  of  a  half-tone  photograph  introduced  in  the  plane  P. . 

Since  we  did  not  have  at  our  disposal  the  elements  to  implement 
exactly  Sayanagui 's  experiment,  we  introduced  in  the  air  space  of  an 
apochroraatic  photographic  lens  a  Dirak  "comb"  consisting  of  an  array  of 
dark  and  clear  parallel  lines  and  provided  the  means  of  rotating  this 
set  of  lines  vri. th  respect  to  the  optical  axis  of  the  instrument. 

Samples  of  photographic  noise  were  manufactured  by  exposing 
Kodak  Kodallth  film  through  a  paper  hand  towel  in  an  ordinary  photo- 
graphic contact  printer.   This  artificial  noise  appeared  to  be  of 
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a  conveniently  random  nature  and  could  be  produced  in  several  densities 
by  varying  the  exposure  time. 

A  periodic  line  pattern  was  superimposed  with  the  noise  to  form 
the  object  and  was  focussed  on  the  image  plane  of  a  lens  containing  the 
filter  described  above.  The  object  and  image  obtained  by  this  procedure 
are  shown  in  Figure  lli. 

The  density  of  lines  tested  as  filters  were  50,  100  and  300  lines 
per  inch.  A  set  of  two  arrays  independently  rotatable  was  also  used. 

The  effect  of  each  filter  was  to  make  the  lens  strongly  aberrant 
astigmatic ally.  Astigmatism  was  maximum  in  the  direction  perpendicular 
to  the  filter  lines  and  zero  in  the  parallel  direction. 

A  single  filter  with  300  lines  per  inch  produced  the  results 
shown  in  Figure  Ik. 

The  good  results  obtained  by  this  method  led  us  to  a  careful 
search  of  the  available  literature  on  the  subject  of  apodization  tech- 
niques for  an  analytical  explanation  of  the  observed  phenomenon  but 
none  seems  to  have  been  published  to  date.  We  believe  that  the  compu- 
tational complexities  of  such  a  theoretical  development  fall  beyond  the 
scope  of  this  research. 

A  final  method  of  filtering  high  spatial  frequencies  of  the  image 
is  by  averaging  the  light  intensity  transmitted  through  the  negative  of 
a  fringe  pattern  in  a  photo-reading  machine.  If  the  scanning  slit  of 
the  instrriment  is  made  large  enough,  the  contribution  of  the  individual 
lines  to  the  overall  intensity  is  diminished  significantly. 

We  have  found  that  to  eliminate  completely  the  influence  of 
the  master  grid  on  the  photo-reading  machine  trace  of  a  fringe  produced 
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id-th  the  master  in  the  image  plane  a  slit  at  least  ten  times  the  size 
of  the  individual  lines  must  be  used. 

For  fringes  produced  with  the  master  in  the  object  plane,  slits 
three  times  wider  than  the  lines  significantly  reduced  their  influence 
on  the  trace.  Figure  15  shows  a  microdensitoraetric  graph  of  the  same 
fringe  scanned  with  a  slit  opening  of  0.081;  mm  for  the  top  trace  and 
0.27  mm  for  the  bottom  one. 

It  is  interesting  to  note  how  all  filtered  fringes,  while  having 
damped  the  image  of  the  grid  lines,  still  reproduce  faithfully  all  varia- 
tions of  the  average  light  intensity.  This  average  is  precisely  what 
carries  the  information  pertaining  to  the  displacement. 


CHAPTER  III 


NUMERICAL  METHODS  FOR  THE  PROCESSING 
OF  miE^   DATA 


3.1  Introduction 

Having  produced  and  recorded  a  moire  fringe  pattern  vdth  the 
photo-optical  instrumentation  analysed  in  Chapter  II,  it  remains  for  us 
to  obtain  the  displacements  and  strains  from  the  continuous  light  inten- 
sity trace  of  the  fringes. 

This  trace,  as  we  said  earlier,  contains  the  useful  signal  given 
by  (2.6)  plus  the  spurious  variations  of  amplitude  and  phase  due  to 
noise. 

We  will  call  the  intensity  trace  function  before  filtering  I  (x) 
so  that 

f'ix)    =   I(x)  +  N(x)  (3.1) 

xdiere  N(x)  is  the  amplitude  and  phase  noise  contained  in  the  signal. 

This  noise  content  will  depend  on  the  amount  of  optical  filter- 
ing performed  on  the  pattern,  however,  for  the  purpose  of  this  chapter, 
no  previous  filtering  will  be  assiimed  to  have  been  made. 

We  are  faced  then  with  the  necessity  of  implementing  a  numerical 
band-pass  filter  to  eliminate  the  unwanted  disturbances  of  the  moire 
light  intensity  trace  given  by  (2.6). 

Once  this  has  been  achieved,  we  desire  to  obtain  the  displace- 
ments given  by  (1.1)  along  the  direction  nornial  to  the  master  grid  lines. 
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The  digitalized  version  of  I  (x)  can  be  obtained  by  the  modifi- 
cation of  a  scanning  microdensitometer  as  the  Joyce  &  Loebl  Mark  III  C 
shovm  in  Figure  16. 

If  this  instrument  is  provided  with  a  potentiometer  arrangement 
that  will  yield  analog  voltages  proportional  to  both  the  density  meas- 
ured on  the  pattern  and  to  the  coordinate  position,  these  voltages  can 
then  be  digitalized  on  a  voltmeter  and  the  numerals  printed  on  tape  for 
convenient  handling  in  a  high  speed  computer. 

The  system,  consisting  of  a  microdensitometer  and  a  digital 
voltmeter,  used  in  the  present  investigation  is  shown  in  Figure  17. 

The  principle  of  operation  of  the  microdensitometer  was  explained 
in  Section  2.7.  It  is  capable  of  performing  linear  density  measurements 
from  0  D  to  6  D,  and  of  linear  magnification  from  1  X  to  1000  X. 

The  digital  data  acquisition  system  is  a  Dymec  1020  B,  that,  at 
the  sampling  rate  of  ten  words  per  minute,  has  a  maximum  resolution  of 
0.0001  V. 

With  this  system,  care  must  be  taken  in  moire  work  to  use  either 
a  non-linear  density  wedge  in  the  microdensitometer,  as  was  described  in 
Section  2.7,  or  to  include  operations  (2.7)  and  (2.8)  in  the  digital 
program  that  will  handle  the  numerical  operations.  Failure  to  do  this 
will  cause  the  non-linear  response  of  the  photographic  film,  used  to 
record  the  fringe  pattern,  to  produce  undesirable  waveform  distortion. 

The  nature  of  the  spurious  disturbances  in  the  moire  light 
intensity  trace  has  been  discussed  in  Section  2.6.  In  order  to  elimin- 
ate these  by  digital  means,  a  successful  band-pass  filter  must  be  imple- 
mented numerically. 
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In  previous  work  [6]  an  in-phase  band-pass  electronic  filter 
was  used  to  eliminate  spurious  variations  in  the  amplitude  and  phase 
of  the  analog  output  of  the  photo-reading  machine.   The  inclusion  of 
the  electronic  filter  produced  transients  of  long  decay  time  that 
lengthened  considerably  the  operating  time  necessary  to  process  the 
fringe  patterns. 

To  overcome  this  difficulty,  numerical  in-phase  and  in-quadrature 
filters  were  combined  in  this  investigation  to  separate  the  amplitude 
and  phase  of  the  light  intensity  trace  and  to  yield  the  phase  informa- 
tion contained  in  the  signal  with  minimum  error. 

The  details  of  these  digital  filters  will  be  presented  in 
Section  3.3. 

Upon  obtaining  the  digitalized  and  filtered  version  of  our 
displacement  curve  (1.1),  we  still  need  to  implement  the  process  of 
derivation  by  which  strains  are  obtained  from  displacements. 

The  definition  of  the  derivative,  as  the  limit  of  a  difference 
quotient,  has  little  value  if  we  possess  only  a  discreet  number  of 
experimentally  determined  values  of  the  function. 

The  ratio  Ay/Ax  becomes  excessively  sensitive  to  small  errors 
in  the  value  of  y  as  Ax  becomes  small  [22].  Therefore,  we  must  resort 
to  a  least  square  approximation  for  a  set  of  values  of  the  function 
in  the  neighborhood  of  the  point  where  we  wish  to  obtain  the  derivative. 

The  influence  of  varying  the  degree  of  the  polynomial  approx- 
imation, and  the  use  of  these  techniques  of  numerical  calculus  for  the 
computation  of  strains,  id.ll  be  seen  in  Section  3.U. 
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If  the  displacement  curve  is  obtained  by  the  use  of  numerical 
in-phase  and  in-quadrature  filters  the  process  of  derivation  can  be 
implemented  in  an  analog  electronic  circuit.  The  details  of  this 
approach  vdll  be  presented  in  Section  3.5. 

To  evaluate  the  merit  of  the  different  procedures  utilized  to 
retrieve  the  phase  information  and  to  derive  the  curve  of  displacements, 
two  auxiliary  tests  were  performed.  These  consisted  of  a  bar  under 
traction  and  a  disk  under  compression.   A  brief  report  on  the  pertinent 
features  of  these  tests  is  given  in  Section  3.2  and  the  comparative 
performance  of  the  different  methods  is  discussed  in  Section  3.6. 

3.2  Auxiliary  Tests 

A  bar  under  axial  tension  and  a  disk  under  diametral  compression 
were  used  as  auxiliary  tests  to  evaluate  the  effectiveness  of  the  numer- 
ical techniques  that  will  be  described  in  this  chapter.  Both  samples 
were  made  of  Hysol  8705  and  engraved  with  a  grid  of  300  lines  per  inch. 

By  means  of  the  tension  test  a  series  of  parallel  equidistant 
moire  fringes  were  produced  that  corresponded  to  the  state  of  homogen- 
eous deformation  in  the  body.  The  ideal,  noiseless,  undistorted  moire 
fringes  would  therefore  give  rise  to  a  pure  cosine  curve  of  constant 
amplitude  and  frequency  as  their  light  intensity  trace.  Our  real 
trace,  after  filtering,  should  approximate  closely  the  ideal. 

The  displacement  curve  along  any  line  parallel  to  the  longi- 
tudinal axis  of  the  sample  should  be  a  line  of  constant  slope  propor- 
tional to  the  strain  in  the  body. 
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For  control  purposes  the  average  strain  in  the  body  was  calcu- 
lated over  a  distance  of  7  in.  along  the  longitudinal  axis  of  the 
sample. 

A  strain  of  3U.85U  x  10~  in. /in.  was  determined.  The  combined 
relative  errors  of  the  microdensitoraeter  location  of  the  first  and  last 
fringe  and  of  the  hand  plot  of  linear  distances  was  0.7%. 

Turning  now  to  the  moire  test  of  the  same  traction  sample,  the 
pattern  was  processed  in  the  system  of  Figure  17.  The  sampling  rate 
was  ten  words  per  second.  The  taped  output  of  the  digital  voltmeter 
was  translated  to  a  punched  card  format  compatible  with  an  IBM  709 
computer. 

The  deck  of  data  cards  was  then  processed  by  the  numerical  and 
analog  methods  described  in  Sections  3.3,   3.ii  and  3.5  to  produce  the 
results  discussed  in  Section  3.6. 

In  a  second  moire  test  a  disk  h   in.  in  diameter  was  loaded 
diametrically,  Figure  18,  and  the  fringe  pattern  was  secured  by  the 
system  of  Figure  17.  The  density  trace  of  the  moire  fringes  across  the 
horizontal  axis  of  the  disk  is  shown  in  Figure  l8. 

In  this  case  a  strong  frequency  modulation  exists  in  the  light 
intensity  variation,  therefore,  a  more  general  state  of  deformation 
than  the  one  existing  in  the  traction  sample  is  obtained. 

The  pattern  was  processed  along  the  x-axis  in  a  manner  entirely 
similar  to  that  described  for  the  traction  sample. 

The  theoretical  solution  to  the  disk  problem  is  well  known  [23] 
and  it  affords  additional  means  of  evaluating  the  accuracy  of  the  numer- 
ical techniques  proposed  here.  The  large  deformations  imposed  here 
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cause  linear  theory  and  the  experiment  to  disagree  somewhat  due  to  the 
relatively  large  deformations  of  the  boundaries  [2k,   25]. 

However,  since  it  was  of  great  interest  to  us  to  determine  if 
the  numerical  procedures  employed  could  retrieve  the  form  of  the  dis- 
placement and  strain  curve  from  the  noise  in  v*iich  they  were  immersed, 
this  test  was  considered  suitable. 

The  density  trace  of  the  fringes  along  the  y-axis  of  the  disk 
is  shown  in  Figure  19.  Under  uniform  illumination  the  amplitude  of  the 
trace  is  seen  to  diminish  as  the  spatial  wavelengths  become  shorter  in 
the  neighborhood  of  the  applied  load. 

Though  the  light  intensity  trace  in  this  direction  was  not  pro- 
cessed to  yield  the  strains,  it  is  shown  here  as  another  good  example 
of  the  spatial  filtering  characteristics  of  the  photo-optical  system 
discussed  in  Chapter  II. 

It  is  clear  that  >dien  the  displacements  along  the  direction 
normal  to  the  master  grid  become  large  enough  to  cause  the  distance 
between  fringes  to  approximate  the  limit  of  resolution  of  the  optical 
system,  the  contrast  in  the  recorded  fringe  pattern  will  be  drastically 
reduced. 

In  the  event  that  this  occurs,  initial  mismatches  of  contrary 
sign  to  the  displacement  expected  in  a  given  region  would  increase  the 
final  spatial  wavelength  of  the  moire  fringes  and  permit  better  analysis. 
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3.}     Displacement  Determination  by  In-Phase 
and  In-Quadrature  Filters 

Having  obtained,  by  the  use  of  the  electro-optical  system,  the 
dlgitalized  version  of  the  moire  light  Intensity  trace  of  the  fringes 
plus  noise,  we  possess  the  numerical  values  of  an  empirical  function. 
We  will  call  this  set  of  numbers  I,  and  the  continuous  function  from 
vrtiich  they  were  obtained  I  (x) . 

We  make  the  following  assumptions  regarding  I  (x) : 

a)  It  is  a  function  that  defines  a  generalized  function 

b)  It  is  band  limited 

c)  The  spectrum  of  the  desired  signal  I(x)  and  that  of 
the  noise  N(x)  are  disjoint. 

The  first  assumption  holds  in  all  cases.  The  other  two  must 
be  taken  as  only  partially  representative  of  the  physical  facts, 

Assiunption  (a)  is  satisfied  by  all  light  intensity  variations 
of  moire  fringes.  In  fact,  I"(x)  satisfied  Dirichlet's  conditions  on 
its  interval  of  definition,  -vrtiich  is  the  length  of  the  model  in  the 
direction  of  scanning. 

Regarding  (b)  the  moire  light  intensity  trace  plus  noise  has 
in  general  a  broad  band  representation  in  the  frequency  plane.  The 
harmonic  content,  however,  can  always  be  cut  short  of  infinity  without 
significant  loss  of  information  contained  in  the  signal. 

Assumption  (c)  may  or  may  not  hold  exactly  depending  on  the 
noise  characteristic  of  the  signal.  If  it  does  not  hold,  the  analysis 
that  follows  will  be  only  approximately  valid,  but  sufficiently  accurate 
for  all  cases  encountered  by  the  writer. 
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If  assumption  (c)  holds,  the  frequency  composition  of  the 
function  to  be  filtered  is  contained  in  a  frequency  interval  f  +  Af/2, 
•where  Af  is  the  band  width  and  f  the  center  frequency.  All  frequencies, 
of  course,  are  spatial  in  moire  work. 

If  we  select  a  sampling  frequency  f  in  our  electro-optical  device 
such  that  f  =  2f  ,  -where   f  is  the  highest  frequency  present  in  the 

S3.  3. 

data,  under  assumptions   (a)    and   (b)  Shannon's  sampling  theorem  permits 
us  to  write 

I,    =    S       a    I.  3.2) 

k       p=  -m    p    k+p 

"■)<■  ^ 

vihere   I  indicates  a  good  approximation  of  I  (x)  obtained  through 


(3.2). 

However,  if  the  coefficients  a  are  the  inverse  Fourier  trans- 
'  P 

*"* 
forms  of  a  particular  transfer  function  in  the  frequency  domain,  I, 

given  by  (3.2)  will  be  a  filtered  version  of  I  (x) .  We  will  call  the 

coefficients  obtained  in  this  manner  a  and  for  the  resulting  filtered 


version  of  I,  we  will  reserve  the  notation  I  ,  thus, 


m 


,-«• 


The  factors  a  become,  in  this  case,  spatial  invariants  that 
play  the  role  of  weight  functions  of  a  linear  filter.  Equation  (3.3)  can 
be  thought  of  as  a  finite  moving  average  process. 

The  basic  relations  of  linear  filter  theory  are  given  in  Appendix 
A  for  convenience  of  the  reader. 

Since  the  information  on  the  displacement  field  is  contained  in 
the  phase  of  the  signal  I(x),  we  wish  to  develop  a  numerical  procedure 
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capable  of  separating  the  phase  variations  from  the  amplitude  changes 
of  a  signal  arriving  at  a  filter. 

A  simple  way  of  accomplishing  this  was  mentioned  by  Gabor  [26] 
and  is  applied  to  our  problem  in  iirfiat  follows. 

Let  us  assume  for  simplicity  that  the  light  intensity  trace  of 
a  moire  pattern  is  a  real  signal  of  the  form 

I(x)  =  a  cos  (Dx  +  b  sin  uux  (3.h) 

This  real  signal  can  be  replaced  by  a  complex  signal  of  the 


form 


Hx)    =  I(x)  +  i  I^(x)  {3.S) 


irtiere  I  (x)  is  defined  as 

q 


hence 


I  (x)  =  a  sin  cux  -  b  cos  uux  (3.6) 


f(x)  =  (a  -  ib)e^'^  (3.7) 


and  then 


1    ,  ,^-1  ^q^-^ 
I(x)  =  Re  $(x)  =  [I(x)^  +  I  (xff    .   e   ^  '^^ 


(3.8) 


q 

The  function  Iq(x)  is  a  signal  in  quadrature  with  I(x)  that 
transforms  the  oscillating  vector  into  a  rotating  vector. 

If  I(x)  is  not  a  simple  harmonic  the  in-quadrature  function  i. 
provided  by  the  Hilbert  transform  of  I(x). 

The  square  brackets  in  (3.8)  represent  the  amplitude  of  I(x) 
vdiile  the  phase  of  I(x)  is  given  by 

1  ^Q^^) 
p(x)  =  tan   .^  (3.9; 
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and  the  amplitude  and  phase  information  contained  in  I(x)  has  been 
separated  to  suit  our  purpose. 

This  transformation  of  an  oscillating  vector  function  in  a 
rotating  vector  function  has  been  vddely  used  in  the  analysis  of  narrow- 
band wave  forms  [27,  28,  29].   The  combination  of  this  transformation 
with  numerical  filtering  has  been  used  by  Ormsby  [30]  and  by  Goodman  [31] 

We  will  refer  in  vrfiat  follows  to  I(x)  as  the  in-phase  signal 
while  I  (x)  will  be  called  the  in-quadrature  signal.   The  latter  is 
obtained  from  the  former  by  a  n/2  phase  shift. 

Since  our  analysis  involves  digital  computation  with  a  signal 

submerged  in  noise,  we  will  be  working  with  the  sampled  version  of 

I  (x)  and  I  (x)  which  we  have  called  I,  and  (I),  . 
q  k      q  k 

The  purpose  of  our  numerical  procedure  will  be  to 

a)  obtain  the  filtered  version  of  I,"  which  we  call  I,  , 

k  k' 

b)  produce  the  90  degrees  out  of  phase  and  filtered 
(I  ),  by  the  expression 

m 

(I  ),  =  E    b  i;  (3.10) 

q  k   p=_ni  P  k+p 

vdiere  the  b  have  the  same  interpretation  as  the  a 
P  P 

in  (3.3),    and  are  derived  from  them  according  to  (3.^:)} 
and  (3.6), 

c)  find  the  guotient  of  (I  )-  and  I  and  introduce  it  in 

,  (i  ), 

0(X)  =  tan"^  -r:^    ,  (3.11) 

^k 
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d)  obtain  the  displacements  given  by 

U(x)  =  po(x),  and  (3.12) 

e)  derive  this  displacement  curve  to  produce  the  strains. 
We  will  apply,  to  fulfill  objectives  (a),  (b)  and  (c),  the 

techniques  proposed  by  Ormsby  and  Goodman,  and  upon  obtaining  a  smooth 
curve  of  phases,  the  change  of  scale  implicit  in  (3.12)  will  yield  the 
displacements. 

The  process  of  derivation  (e)  will  be  considered  in  Section  3.h. 

(!)     Goodman's  linear  combination  of  Tukey  filters. 
The  ideal  transfer  function  of  an  in-phase  and  in-quadrature 
filter  is  sl^own  in  Figure  20. 

Goodman  [31]  proposed  a  linear  combination  of  Tukey  filters 

that  is  centered  at  f   is  essentially  flat  from  f  -  n/m  to  f  +  n/m 

°  0  0       ' 

and  zero  outside  the  interval  f  -  3n/m  to  f  +  3v/m.      The  coefficient 
m  is  the  number  of  terms  in  the  trigonometric  series  of  the  Tukey 
filter  transfer  function. 

The  filter  function  that  results  from  the  linear  combination  is 
shown  schematically  in  Figure  21. 

The  filtering  finite  moving  in-phase  and  in-quadrature  average 
processes  in  the  space  domain  are  given  by 


^k 


)    o  m-1  (0.5U  +  0.16  cos  j  -)  sin  i  -  cos  i  -  cos  i  f 

ii  +  £  Y     —— _JI! __£! __i" "'  0 

"   "  j=l  J 

U     8  """"^ 
^m  \  "  m  J?!  (°-^^  *  0-^6  cos  j  ^)  cos  j  ^  cos  j  i:^{l^^^^l^_.) 


+  -   (O.5I4  +  0.ii6  cos  tt)  cos  tt  cos  m  f  (I    +1   )] 
"^  o  k+m    k-m 


(3.13) 
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'V. 


,    ra-1   (0.51;  +  0.i|6  cos  J  -)    sin  2j  -  sln^  j    f 


ii     E     


.8  '"-^ 


-1 


^m    -^n    ^°-^^  ■*■  °-^^  ^°^  ^  ^)    ^°s  j  -  sin  J   f  (I,     .   +  I        ) 

+  -  (O.5I4  +  0.ii6  cos  n)    cos  n  sin  m  f   (I,         +1        )]. 
f"  0     k+m         k-m   -^ 

(3.1a) 

The  values  of  the  filter  transfer  function  in  the  interval 
(f^  +  n/m|  to  |f^  -  n/m|  is  called  the  filter  roll-off.  For  Goodman's 
filter  the  roff-off  is  a  straight  line  inclined  vdth  respect  to  the 
vertical  axis.  In  this  interval  the  real  and  ideal  filter  differ  appre- 
ciably. The  real  filter  allows  noise  of  frequency  close  to  that  of  the 
signal  to  pass  onto  the  output. 

Real  filters  vary  as  to  the  monotonlc  function  that  is  selected 
to  connect  the  horizontal  portion  of  the  transfer  function  with  the 
frequency  axis  [32].  Tukey,  as  we  said  above,  selected  a  straight  line. 

This  produces  a  sharp  discontinuity  of  the  filter  function  at 
^o  ±  3TT/m  which  in  turn  causes  any  series  approximation  used  to  express 
analytically  the  function  to  overshoot  the  value  one,  due  to  Gibbs' 
effect,  at  the  discontinuity. 

Another  undesirable  feature  of  this  filter  transfer  function 
is  the  dependence  of  the  band  width  Af,  the  roll-off  frequency  interval, 
and  the  number  of  terms  in  the  finite  moving  average  process  on  a  single 
parameter  m. 
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The  selection  of  one  coefficient  to  suit  these  three  different 
purposes  restricts  the  applicability  of  this  particular  filter  to  prob- 
lems where  the  signal  has  a  narrow-band  frequency  composition  and  the 
noise  and  signal  spectrum  are  disjoint. 

The  only  other  independent  parameter  of  the  filter  is  the  cen- 
tral frequency  f  . 

The  practical  implications  of  the  comments  made  above  will  be 
explored  further  in  Section  3.5  where  comparative  quantitative  results 
are  presented. 

(2)  Onnsby's  filter. 

A  schematic  representation  of  an  Ormsby  band-pass  filter  trans- 
fer function  is  shown  in  Figure  22. 

The  roff-off  function  of  this  filter  is  composed  of  two  para- 
bolic shoulders  at  f  and  f  connected  by  a  straight  line.  The  length 

C        X" 

of  the  interval  of  definition  of  the  parabolic  transition  can  be  modi- 
fied by  varying  the  value  of  the  parameter  k. 

The  independent  parameters  of  the  filter  are  the  number  of 
terms  in  the  trigonometric  series  approximation  of  the  transfer  function 
N,  the  center  frequency  f  ,  the  cut-off  frequency  f  ,   the  largest  fre- 
quency f. ,  and  the  value  of  k  explained  above  which  is  subject  to  the 

condition 

f ,  -  f 
k  <  ^  g  ^  (3.15) 

The  parabolic  transition  increases  the  rate  of  convergence 
of  the  trigonometric  series  to  the  transfer  function  in  the  neighbor- 
hood of  f  ,  thus  permitting  better  filtering  of  noise  close  to  the 


hh 


signal  spectrum.  The  independence  of  f .  and  f  insures  a  narrow  roll- 

x*  c 

off  frequency  interval. 

The  transfer  function  of  this  filter  is 


H(f)  =  ^ 
1 


if 


if 


1  -  a(f  -  f  )  if 
^            c 

1  -  a(f  -  f     .2  if 

bf  +  c  if 

bf  +  c  if 


a(f  -  f^)        if 
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f        <    f    <    f     +  kAf 
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f       +  kAf,     <    f     <    f 
C]_  1'         '    '  ci' 


f     +  kAf       <    f    <     f^   -  kAf, 
c  1'        '  '    t  1' 


f^^  +  kAf^l    <  lf|    <    |f^^  -  kAf^l 

f^     -  kAf^l    <  |f|    <   |f^| 

f.J  <  kl  <  If.   +  kAf  I 
ti  ti  1 


(3.16) 


where  the  parameters  a,  b,  and  c  are  solved  in  terms  of  k  to  satisfy 
the  condition  (3.16). 

For  the  above  transfer  function  the  in-phase  and  in-quadrature 
finite  moving  average  processes  that  yield  the  filtered  versions  of 

I  and  (I  )  can  be  found  to  be 

m       q  m 
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E    h  cos  2TTn  X     I 


(3.17) 


N 

(I  )  =  E   h  sin  2TTn  X  l"" 
m  q   n=  -N  "         o  n+m 


(3.18) 
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where 

sin  2TTn[\  +(l-k)X   ]+  sin  2TTn(X  +kX   )-sin  2TTn(X  +X   )-sin  2TTn  X 
h     -  "•   c  r  c       r  c     r  c 


n 


Uk(l-k)X^(TTn)^ 


(3.19) 
and 

K-k'         't'S'    '="ft'    'r-\-\  (3.20) 

The  filters  (3.13),  (3.1Ii),  (3.17)  and  (3.18)  were  implemented 
in  an  JEM  709  programming  code.  By  several  trial  runs  -with  different  . 
value  of  the  independent  parameters  the  optimum  performances  of  both 
filters,  acting  on  the  data  from  the  two  auxiliary  tests  described  in 
Section  3.2,  were  obtained.  The  quantitative  comparison  of  their  per- 
formances is  presented  in  Section  3.S. 

After  optimization  (3.11)  and  (3.12)  were  implemented  to  yield 
the  displacements.  The  corresponding  programs  are  included  in 
Appendix  C. 

3.U  strain  Determination  by  Numerical  Differentiation 
>     ~~~~~~""   — — — — ^-^— ^— — ^— — — ^^— ^^— — — — 

When  the  displacement  curve  has  been  obtained  in  the  manner 
described  above,  it  remains  for  us  to  derive  numerically  this  empirical 
function  to  produce  the  strains. 

Previous  to  this  research  [37],  the  use  of  an  in-phase  electronic 
band-pass  filter  did  not  produce  a  displacement  curve  sufficiently  smooth 
for  derivation  purposes  and  a  local  averaging  process 

\  =  7  ^   6  \  (3.21) 

i=k-^ 


he 


had  to  be  applied  twice  before  differentiation.  The  summation  extended 
over  the  number  of  data  points  in  one  fringe. 

The  resulting  strains  were  calculated  by  a  finite  difference 
approximation 

,    \^1   -  \-l 

\   =  ^ (3.22) 

k 

where  A  was  the  sampling  wavelength. 

The  strains  obtained  in  this  manner  were  again  averaged  twice 

by  the  formula 

\   =  7  .^,  6  \  (3.23) 

The  analysis  of  this  procedure  reveals  several  iindesirable 
features  which  we  will  now  discuss. 

The  necessity  of  an  amplitude  normalization  has  already  been 
mentioned  as  a  drawback  of  the  method.  Furthermore,  since  no  weights 
are  associated  with  the  point  values  of  the  function  in  the  neighbor- 
hood of  where  the  derivative  is  sought,  and  the  finite  difference 
approximation  is  sensitive  to  small  errors  in  these  values,  the  con- 
vergence of  the  averaging  process  to  the  mean  is  slow. 

Though  repeated  application  of  (3.21)  and  (3.23)  does  smooth 
out  the  displacement  and  strain  curve,  it  causes  us  to  lose  information 
at  the  beginning  and  end  of  the  moire  light  intensity  trace.  In  fact, 
each  average  shifts  the  first  point  at  which  the  smoothed  curve  can  be 
defined  a  distance  6/2  along  the  direction  of  scanning. 
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A  series  of  two  averaging  processes  for  the  displacements  and 
two  for  the  strains  would  erase  all  information  on  the  strain  curve 
contained  in  the  two  moire  fringes  closest  to  every  boundary. 

In  the  thermal  tests  of  Chapter  IV,  the  strain  determination  at 
the  boundaries  of  the  model  was  critical  since  the  largest  circ\imfer- 
ential  stresses  occurred  at  these  points.  To  overcome  this  difficulty 
and  to  obtain  increased  accuracy  in  the  derivation  procedure,  we  era- 
ployed  a  method  suggested  by  Lanczos  [22].  It  provides  for  the  approx- 
imation, in  the  least  square  sense,  of  a  polynomial  of  arbitrary  degree 
n,  to  a  set  of  experimentally  obtained  values  of  a  function  in  the 
neighborhood  of  the  point  where  we  wish  to  define  the  derivative.  This 
polynomial  is  then  derived. 

The  expression  for  this  derivative  is 


s 

E   Qt)  (x+Q-h) 

dp(x)  ^  g=  -s 

dx        ^       0 
2  E  cc  h 


(3.2U) 


^=1 

where  h  is  the  sampling  wavelength  of  the  data  and  s  the  degree  of 
the  polynomial. 

Since  the  least  square  approximation  is  performed  in  the  small, 
no  significant  error  is  committed.  The  variation  of  the  degree  of  the 
polynomial  can  be  conveniently  performed  by  the  inclusion  of  s  as  an 
independent  parameter  in  the  programming  code. 

To  obtain  the  derivative  at  the  beginning  and  end  of  the  dis- 
placement curve,  (3.2li)  can  be  modified  [22],  or  a  reflection  of  the 
experimental  curve  can  be  performed  with  respect  to  the  first  and  last 
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data  point.  Both  methods  have  been  used  in  the  numerical  examples 
of  Section  3.5. 

The  fitting  of  a  polynomial  to  the  entire  displacement  point 
function  by  the  method  of  least  squares  did  not  offer  a  satisfactory 
solution  because  it  implied  making  a  decision  on  the  degree  of  the  poly- 
nomial approximation  that  endangered  the  objectivity  of  the  procedure. 

The  techniques  of  differentiation  that  rely  on  integration  were 
considered  [22,  3},   }k] .     The  one  proposed  by  Lanczos,  as  the  limit  of 
{3.2k)  when   the  sampling  wavelength  becomes  very  short,  yielded  the 
best  results  with  less  computing  time. 

At  the  limit  (3.2ii)  becomes 

^--hl       ^   ^(^■^^)  dt  (3.25) 

•Mhere   dx  is  the  sampling  wavelength  and  T]  defines  the  limits  of  the 
integration  interval. 

The  strain  curve  obtained  by  (3.2J4)  or  (3.25)  may  or  may  not 
be  smooth  enough  for  our  purposes  depending  on  the  amount  and  the  kind 
of  noise  present  in  our  signal.  Since  the  derivative  process  tends  to 
magnify  small  errors  in  the  displacement  function  that  is  being  derived, 
we  may  wish  to  apply  a  final  smoothing  process  to  the  strain  curve 
obtained.  For  this  purpose  we  select  a  filter  of  the  form  [3$] 

i 

e '  =  ^  .  w  e(n+r)  (3.26) 

"   r=  -X  ^ 

vdiere  e^  are  the  raw  values  of  the  derivative,  w  the  weights  of  the 


r 


filter  and  e  the  smooth  derivative  values. 
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If  w  =  •^r  for  all  1*5  (3.26)  reduces  to  a  local  average  process 
identical  to  (3.23),  however,  (3.26)  affords  the  possibility  of  select- 
ing the  w  ' s  to  produce  a  more  rapid  convergence  of  the  noise  ridden 
function  to  the  smooth  for  smaller  values  of  i.  The  choice  of  the  best 
w  for  a  given  filtering  application  depends  on  the  type  of  noise  in 
the  curve  of  derivative.  As  vdll  be  seen  in  Section  3.6,  several  trials 
are  often  necessary  to  obtain  the  optimum  filter  weights. 

The  programming  codes  corresponding  to  this  section  are  pre- 
sented in  Appendix  C . 

3.5  Strain  Determination  by  a  Hybrid  Computer  System 

The  previous  two  sections  were  concerned  with  purely  numerical 
methods  of  displacement  and  strain  determination.  We  will  now  describe 
a  process  that  combines  the  digital  procedure  of  in-phase  and  in- 
quadrature  filtering  that  produces  the  displacement  curve  with  a  filter- 
ing and  derivative  circuit  obtained  on  an  analog  computer.  For  this 
purpose  we  combine  a  Ramo  Wooldridge  RW  300  digital  computer  with  an 
Applied  Dynamics  AD2  digital  system  to  produce  a  hybrid  computer. 

Operations  (3.3),  (3.10),  (3.11)  and  (3.12)  are  implemented  and 
performed  on  the  digital  section  of  the  hybrid  vftiile  the  equivalent  of 
Equations  (3.2U)  or  (3.25)  are  done  on  the  analog  section  of  the  system. 

The  displacement  curve  produced  by  digital  methods  is  smoothed 
by  an  analog  exponential  filter  circuit.  This  filtered  version  of  the 
displacement  curve  is  derived  by  another  circuit  to  produce  the  strain 
distribution.  The  details  of  these  circuits  are  given  in  Section  3.6. 
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The  advantage  of  this  alternate  process  of  differentiation  is 
the  speed  with  which  the  curve  of  strains  can  be  obtained  from  the 
analog  and  the  ease  with  which  the  parameters  that  control  the  amount 
of  filtering  on  the  displacement  and  strain  curve  can  be  varied. 

We  see  that,  though  the  digital  methods  of  Sections  3.3   and  3.1| 
permit  us  to  take  advantage  of  the  continuous  law  (1.1),  reduce  the 
real  time  of  pattern  analysis  and  increase  the  objectivity  of  moire 
data  processing,  the  versatility  of  the  procedure  is  somevrfiat  less  than 
ideal  since  several  runs  must  be  made  to  optimize  the  independent  para- 
meters involved.  This  can  make  the  method  costly  in  computer  time. 

As  far  as  we  have  been  able  to  determine,  there  exists  to  date 
no  electronic  system  that  can  perform  the  orthogonalizations  (3.6)  or 
(3.10),  so  we  have  had  to  settle  for  the  combination  of  numerical  and 
analog  data  processing  described  above. 

3.6  Comparative  Performance  of  the  Proposed  Techniques 

In  this  section  we  will  present  the  results  obtained  by  the  use 
of  the  numerical  and  analog  techniques  discussed  in  Sections  3.3,   3.h 
and  3.5  acting  on  the  digitalized  version  of  the  light  intensity  trace 
corresponding  to  the  tests  of  Section  3.2, 

For  the  sake  of  brevity  we  will  agree  to  call  this  section 
Filter  1  to  Goodman's  linear  combination  of  Tiikey  filters,  and  Filter  2 
to  Ormsby's  filter. 

a)  Displacement  determination 

As  was  said  in  Section  3.3,  (3.13)  and  {3.1k)   that  correspond 
to  Filter  1  were  implemented  in  an  IBM  709  programming  code. 
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The  independent  parameters  of  this  filter  are  the  number  of  terms  in 
the  trigonometric  approximation  of  the  transfer  function  m,  and  the 
central  frequency  of  the  signal  f  . 

In  another  program  (3.17),  (3.18)  and  (3.19)  were  implemented 
for  Filter  2.  This  filter  has  as  independent  parameters  the  number 
of  terms  in  the  trigonometric  expansion  N,  the  center  frequency  f  ,  the 
cut-off  frequency  f  ,  the  termination  frequency  f .  and  k  a  factor  that 

C  u 

defines  the  extent  of  the  parabolic  transition  of  the  roll-off.  Both 
programs  contained  (3.11)  and  (3.12)  to  produce  the  displacements. 

Traction  test  results 

The  homogeneous  state  of  deformation  of  the  model  affords  the 
possibility  of  determining  quite  accurately  the  band  width  of  the 
numerical  filter.  In  fact,  the  transfer  function  of  the  filter  reduces 

to 

1  for  f  =  f 


( 
H(f)  ={ 


0  for  f  /  f 

V  0 

However,  the  actual  intensity  trace  of  the  moire  fringes  of  the  trac- 
tion sample  showed  that  at  a  constant  sampling  speed  the  number  of 
data  points  per  fringe  varied  somev±iat  due  to  noise.   Therefore,  to 
determine  the  parameter  f  the  average  value  of  the  wavelengths  was 

chosen.  This  f  was  the  same  for  both  Filter  1  and  Filter  2. 
o 

The  criterion  for  the  selection  of  the  number  of  terms  in  the 
trigonometric  expansion,  m  and  N  for  Filters  1  and  2,  is  different  for 
each  filter. 
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Since  for  Filter  1,  in  defines  also  the  band  width  of  the 
transfer  function,  it  is  necessary  for  it  to  be  large  enough  so  that 
the  quotient  Tr/m  be  very  small.  This  increases  the  number  of  terms  in 
(3.13)  and  (3.1ii),  making  the  use  of  the  filter  costly  in  computer  time' 
to  retrieve  constant  phase  cosine  waves  immersed  in  noise. 

For  Filter  2  the  situation  is  entirely  different.  The  exist- 
ence of  independent  parameters  to  define  the  number  of  terms  in  (3.17) 
and  (3.18)  and  the  extent  of  the  band  width  make  it  very  practical  for 
the  purposes  of  this  test. 

The  roll-off  frequency  interval  can  also  be  narrowed  conven- 
iently by  the  use  of  the  parameters  f  ,   f .  and  k. 

The  displacement  curve  obtained  with  Filter  2  is  shown  in 
Figure  23.  The  computer  values  appear  as  discrete  points  -id^ile  the 
theoretical  curve  is  drawn  in  full  line. 

The  displacements  were  obtained  over  base  lengths  of  measure- 
ment of  6,  6/1;,  6/5,  6/8,  6/10  and  6/iiOO.  The  standard  deviations  at 
the  different  levels  of  fractional  fringe  interpolation  varied  between 
an  upper  bound  of  +  0.001  p  or  21'  36"  of  arc  and  a  lower  bound  of 
+  0.0003  p  or  6'  29"  of  arc.  The  maximum  variance  between  two  means 
at  different  levels  of  interpolation  was  found  to  be  0.00025  p. 

The  consistency  of  these  results,  and  the  fact  that  no  cyclic 
variation  in  the  errors  of  the  displacement  detennination  present  in 
earlier  work  [6]  was  found  to  exist  in  this  method,  demonstrate  that 
a  substantial  improvement  has  been  achieved. 

The  results  reported  above  indicate  that  the  minimum  instantan- 
eous angle  that  can  be  meaningfully  measured  by  this  procedure  is  1  . 
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This  implies  that  we  have  effectively  divided  the  grid  pitch  into  360 
equal  parts,  that  is,  a  five-fold  increase  over  previous  methods  of 
fractional  moire  fringe  interpolation. 

Disk  under  compression  test  results 

In  this  case  the  light  intensity  trace  along  the  x-axis  of  the 

ring  presents  strong  spatial  frequency  modulation,  Figure  18.  This 

indicates  that  the  frequency  composition  of  the  signal  is  contained  in 

an  interval  f  +  Af/2,  where  Af  is  not  negligibly  small  as  in  the  case 
o 

of  the  traction  sample.  This  fact  allows  Filter  1  to  perform  more 
effectively  here  than  in  the  previous  test  and  we  will  report  the  com- 
parative performance  of  both  filters. 

The  central  frequency  f  was  detennined  in  terms  of  data  point 
numbers  by  inspection  of  a  standard  IM  80  x  80  listing  of  the  data. 
For  Filter  1  several  trial  runs  had  to  be  made  to  obtain  the  best  combi- 
nation of  the  only  two  parameters  available,  m  and  f  . 

Turning  now  to  Filter  2,  it  again  proved  to  be  easier  to  tune 
in  this  test  than  Filter  1.  The  reasons  for  this  were  explained 
earlier.  The  displacement  curve  obtained  is  shown  in  Figure  2li.  Both 
filters  proved  to  be  able  to  retrieve  the  form  of  the  displacement 
curve  precisely.  The  good  approximation  between  the  numerical  values 
of  the  strains  derived  below  and  those  obtained  by  independent  means, 
confirms  quantitatively  this  accuracy. 
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b)  Strain  determination 

To  compute  the  strains  from  the  displacement  curve  produced  by 
the  filters,  (3.2U),  (3.25)  and  (3.26)  were  programmed  and  the  compara- 
tive performance  of  {}.2k)   and  (3.25)  analysed.  Due  to  the  wealth  of 
data  points  in  our  possession,  (3.25)  could  be  used  with  some  advantage 
over  (3.2U). 

To  obtain  the  derivatives  at  the  beginning  and  end  of  the  sample 
it  was  found  useful  to  reflect  the  displacement  curve  of  the  disk  \inder 
compression  with  respect  to  the  first  and  last  data  point  locations. 
However,  for  the  displacement  curve  of  the  traction  sample  the  reflec- 
tion technique  produced  a  sharp  discontinuity  with  no  gain  in  the  accu- 
racy of  the  derivative.  In  this  case  the  weights  of  the  polynomial 
approximation  implicit  in  (3.2li)  were  varied  at  the  ends  of  the  trace 
as  recommended  by  Lanczos  [22]. 

For  the  strain  curve  obtained  from  the  traction  test  no  smooth- 
ing  with  (3.25)  was  performed,  so  the  values  reported  below  correspond 
to  the  straightforward  application  of  (3.2ii). 

For  the  ring  under  compression  test  (3.25)  was  applied  since 
some  scatter  of  values  was  noticeable  after  derivation. 

Traction  test 

The  strains  obtained  in  the  manner  described  above  between  two 
successive  fringes  are  shown  in  Figure  25a.  The  small  scale  of  the 
plot  is  intended  to  indicate  the  variations  about  the  mean  and  the 
average  strain  determined  in  Section  3.2. 
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A  histogram  was  prepared  to  show  the  relative  frequency  of 
occurrence  of  the  random  variable.  This  variable  was  subdivided  into 
discrete  intervals  spanning  10  x  10"  in. /in.   The  resulting  graph  is 
shown  in  Figure  25b. 

The  characteristics  of  this  distribution  are  of  interest. 
A  saddle-shaped  function  replaces  the  more  familiar  bell- shaped 
Gausian  distribution  we  might  have  expected  to  obtain.  However,  the 
saddle  shape  insures  that  98^  of  the  values  of  the  random  variable  fall 
within  the  interval  u  +  50  x  10"  vrtiich  is  quite  sufficient  for  most 
engineering  applications. 

The  individual  strains  plotted  in  Figure  25a  and  included  in  the 
listing,  correspond  to  a  base  length  of  measurement  of  6/1^00.  If  the 
3%%   certainty  that  any  individual  measurement  is  contained  in  the  inter- 
val M,  +  50  X  10   in. /in.  is  taken  as  a  figure  of  merit  of  the  proposed 
techniques,  we  obtain  an  accuracy  of  SO   x  10"  in. /in.  in  the  strain 
determination.  This  is  well  in  excess  of  previous  results  in  moire  work. 

Disk  under  compression  test 

The  numerical  differentiation  of  the  displacement  curve  of 
Filters  1  and  2  were  performed  with  (3.2lt)  and  subsequently  filtered 
with  (3.26).  Several  trial  runs  were  necessary  to  obtain  the  best 
values  of  the  parameters  s,  i  and  w  .  The  resulting  curve  of  strains 
can  be  seen  in  Figure  26,  where  one-half  of  the  strain  curve  corre- 
sponds to  Filter  1  and  the  other  to  Filter  2,  for  the  same  values  of 

s,  i  and  w  . 
'       r 
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The  full  line  curve  corresponds  to  the  moire  hand  analysis  of 
the  problem  done  by  Sciaramarella  [1],  while  the  dashed  curve  represents 
the  results  obtained  by  Durelli  and  Mulzet  [2U]. 

The  latter  investigators  reported  natural  strains  defined  as 

e   =  ^n(l  +  \      °) 

Tijhere  i     and  ^„  are  the  initial  and  final  lengths,  respectively, 
o      1 

Since  moire  yields  Eulerian  strains  defined  as 

the  results  of  [2ij]  were  transformed  accordingly  to  permit  comparison. 

The  value  of  Poisson's  ratio  for  both  [l]  and  our  tests  was 
0.U65  while  in  [2U]  this  value  was  between  O.U?  and  0.U8.  Since  the 
strain  distribution  is  a  function  of  this  ratio  some  small  inaccuracy 
has  occurred  [2U]. 

The  modulus  of  elasticity  in  [1]  and  here  is  Eulerian  while  in 
[2li]  it  is  natural;  however,  for  strains  less  than  0.03  in. /in.  both 
almost  coincide  [2U]. 

The  percentage  of  errors  committed  by  Filters  1  and  2  with  respect 
to  the  strains  obtained  in  [1]  and  [2li]  is  presented  in  the  following 
table. 

As  can  be  seen  from  the  table  and  Figure  26,  Filter  2  produces 
a  smoother  curve  that  more  nearly  approximates  those  of  [1]  and  [2l|]. 
Both  filters  do  not  perform  well  near  the  edges  of  the  disk  where  the 
strains  are  approaching  zero.  This  is  believed  to  be  due  to  the  fact 
that  long  wavelengths  have  not  been  properly  filtered  in  the  quadrature 
process. 
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TABLE  1 

RELATIVE  ERRORS  IN  THE  STRAINS  OF  THE  DISK 
UNDER  COMFRESSION 


Filter  1 

Filter  2 

ication 

Relative 

to 

Location 

Relative 

to 

x/r 

[1] 

[21;] 

x/r 

[1] 

[2U] 

% 

% 

% 

% 

-1.0 

0 

0 

0.0 

i;.3 

s.^ 

-0.9 

W4.5 

UU.5 

0.1 

0.0 

2.7 

-0.8 

20.8 

20.8 

0.2 

1.6 

1.1 

-0.7 

22.0 

17.0 

0.3 

7.5 

0.0 

-0.6 

0.0 

0.0 

o.U 

9.3 

0.0 

-0.5 

10.7 

ii.O 

0.5 

5.8 

3.U 

-O.k 

9.h 

0.1 

0.6 

16.0 

16.0 

-0.3 

16. U 

9.1 

0.7 

11.0 

16.0 

-0.2 

3.7 

0.0 

0.8 

21.0 

21.0 

-0.1 

1.1 

3.8 

0.9 

60.0 

60.0 

0.0 

3.2 

I4.2 

1.0 

0.0 

0.0 

Now  we  wish  to  present  the  results  obtained  with  the  hybrid 
system  of  Section  3.5. 

The  phases  obtained  with  Filter  2  were  smoothed  with  the  cir- 
cuit of  Figure  27a;  the  resulting  displacement  curve  is  shown  in  Fig- 
ure 27b.  These  displacements  were  derived  by  the  circuit  of  Figure 
27c  to  yield  the  strains  of  Figure  27d.  The  disymmetry  in  the  strain 
curve  is  due  to  the  fact  that  the  digital  and  analog  sections  of  the 
hybrid  system  are  time  locked;  thus,  there  is  a  progressive  delay  in 
the  rate  at  which  the  data  stored  in  the  RW  300  is  relayed  to  the  AD2. 

The  resulting  strain  curve  compares  with  that  obtained  in  [l] 
as  shown  in  Figure  28.  Again  the  maximum  strain  at  the  center  of  the 
ring  is  3.6^  lower  than  that  of  [1].  Several  sections  of  the  curve 
coincide  to  within  +  0.5^.   Other  errors  calculated  are  %.3%   at 
x/r  =  0.3,  and  8.6^  at  x/r  =  0.50. 
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From  the  consideration  of  Figures  26  and  28  we  conclude  that 
the  hybrid  system  yields  a  strain  curve  with  less  error  than  that  pro- 
duced by  purely  numerical  techniques  for  the  amount  of  filtering  used 
in  the  comparison.  However,  further  numerical  filtering  could  con- 
ceivably diminish  the  errors. 

The  strains  produced  by  analog  means  also  present  larger  errors 
at  the  edges  of  the  disk  for  the  same  reasons  explained  above. 

J-ft  might  appear  paradoxical  that  at  this  final  stage  of  the 
development  of  a  completely  computerized  system  of  strain  analysis 
we  should  check  our  results  against  those  produced  by  hand  analysis. 
However,  it  must  be  remembered  that  the  ability  of  the  human  eye  to 
detect  trends  in  a  curve  submerged  in  noise  and  to  select  the  proper 
tangent  is  often  difficult  to  replace  by  another  system. 

In  a  case  such  as  this  of  a  ring  under  compression  where  an 
intelligent  analyst  has  a  good  idea  of  what  results  to  expect,  the 
fact  that  any  piece  of  myopic  hardware  can  compete  idth  him  is  quite 
gratifying  since  total  objectivity  and  tremendous  savings  in  labor  and 
time  are  welcome  bonuses. 


CHAPTER  IV 


THE  THERMAL  TEST 


i;.l  The  Experimental  Set -Up 

A  general  view  of  the  experimental  set-up  can  be  seen  in 
Figure  29.  Basically  it  consists  of  four  units:   the  model,  the  photo- 
graphic equipment,  the  heating  and  cooling  system  and  the  temperature 
control. 

a)  The  models 

The  models  were  disks  made  of  AISI  30I4  stainless  steel  of 
nominal  composition  0.06^  C,  16^  Cr,  and  8^  Ni.  The  disks  were  stress 
relieved  by  heat  treating.  The  outer  diameter  of  all  models  were 
5.5  in.,  the  interior  diameter  varied  from  0.5  to  2.5  in.  in  steps  of 
0.5  In.  The  thickness  of  all  rings  was  1  in. 

We  will  report  here  only  the  results  corresponding  to  what  we 
will  call  Ring  No.  1  with  OD  =  5.5  in.  at  1120°F  and  ID  =  2.S   in.  at 
920°F,  and  Ring  No.  2  with  OD  =  S.S   in.  at  l58o'¥  and  ID  =  0.5  in.  at 
900°F. 

The  models  were  engraved  with  a  chemically  etched  cross-grid 
of  300  lines  per  inch  on  one  face.  A  description  of  the  engravure 
technique,  together  with  other  information  regarding  the  preparation 
of  models  for  high  temperature  testing  can  be  found  in  Appendix  A. 
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b)  The  photographic  equipment 

The  photographic  equipment  consisted  of  a  Saltzman  enlarger 
converted  to  a  photographic  camera.  The  details  of  the  transformation 
were  reported  earlier  by  Ross  [36].  The  interesting  features  of  this 
camera  are  its  rigidity,  the  precision  gearing  that  permits  accurate 
positioning  of  its  back  and  front,  and  the  microinch  counters  that  allow 
faithful  reproduction  of  relative  positions. 

Perfect  parallelism  of  the  image  and  object  planes  can  be 
obtained  and  preserved  for  a  vide  range  of  object  and  image  distances. 
The  back  of  the  camera  can  be  rotated  360  . 

The  film  was  placed  against  the  emulsion  side  of  a  high  resolu- 
tion photographic  plate  vrtiich  served  as  the  master.  Through  this  plate 
the  film  was  exposed.  The  moire  fringe  pattern  was  formed  therefore  in 
the  image  plane  of  the  camera  between  the  image  of  the  model's  lines 
and  the  master  plate. 

To  make  adjustments  of  focus  and  to  obtain  the  desired  mismatches, 
the  image  on  the  focal  plane  of  the  camera  could  be  inspected  with  the 
aid  of  a  removable  ground  glass.  Once  the  adjustments  were  made  the 
ground  glass  would  be  removed,  the  film  loaded  in  its  place,  and  the 
camera  closed. 

The  convenient  mismatches  could  be  obtained  by  simply  changing 
the  distance  between  the  front  and  back  of  the  camera,  while  the  object 
was  kept  focussed.  The  rotation  of  the  back  of  the  camers  permitted  the 
elimination  of  all  rotation  mismatch  by  minimizing  the  number  of  fringes 
in  the  image.  A  further  verification  of  the  absence  of  rotation  mismatch 
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could  be  performed  by  varying  the  size  of  the  image.  If  the  number  of 
fringes  increased  or  diminished  without  changing  direction  no  such  mis- 
match was  present. 

The  apochromatic  photographic  lens  with  a  time  shutter  tested 
in  Chapter  II  was  used  to  form  the  image. 

A  variety  of  films  were  evaluated  prior  to  running  the  test. 
The  films  found  to  record  better  the  fringes  at  different  temperatures 
are  given  below. 

TABLE  2 
Fim  PERFORMANCE  AS  A  FUNCTION  OF  TEMPERATURE 


Maximum 
Temperature 
(°F) 

Film 

ASA 
Speed 

Developer 

200 

Super  Panchro 
Press 

250 

DK  50 

300 

Royal  Pan 

1400 

DK  50 

500 

Contrast 
Process  Ortho 

so 

DK  11 

800 

Super  Speed 
Ortho 

1^5 

DK  50 

1600 

Royal  Ortho 

Uoo 

DK  50 

As  we  can  see,  for  low  temperatures,  panchromatic  films  record 
well  the  moire  patterns.  As  temperatures  increase,  the  fringe  contrast 
is  lost  with  this  type  of  films  and  orthochromatic  films  should  be  used. 
The  reason*  for  this  v/as  given  in  [37]. 
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An  important  factor  in  the  successful  photography  of  fringes  at 
high  temperatures  has  been  found  to  be  the  type  of  light  sources  used 
[37].  For  all  our  photography  an  actinic  cold  light  source  with  pre- 
dominant radiation  in  the  wavelength  X   =  0.5  x  10~  cm  was  used. 

Good  contrast  was  obtained  in  the  temperature  range  between 
300  F  and  1600*T''  with  several  orthochromatic  films  of  varying  speeds. 
However,  the  fastest  should  be  preferred  due  to  the  convectlve  air 
currents  that  may  exist  in  the  space  between  the  model  and  the  lens. 
These  changes  in  the  air  density  caused  by  the  heat  gradient  may  seri- 
ously hinder  photography  by  producing  what  appears  to  be  a  movement  of 
the  object. 

To  keep  them  at  a  minimum,  the  model  was  surrounded  by  a  chamber 
and  evacuated  by  a  vacuum  pump.  A  polished  plate  of  fused  silica 
separated  about  1/U  in.  from  the  model  surface  closed  the  chamber  towards 
the  camera  providing  an  excellent  window.  Finally,  about  $   in.  above 
the  model  surface,  a  polished  glass  plate  was  placed  parallel  to  the 
fused  silica. 

With  these  precautions  and  fast  photographic  films  the  convec- 
tion current  problem  was  overcome  completely. 

c)  The  heating  and  cooling  system 

Th^e  induction  heating  technique  was  selected  because  it  allowed 
an  accurate  control  of  the  heat  input  and  a  perfect  simulation  of  the 
boundary  conditions  of  the  problem. 

The  induction  heating  coils  were  wrapped  around  the  outer  circum- 
ference of  the  rings.  The  diameter  of  the  coils  exceeded  that  of  the 
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rings.  The  diameter  of  the  colls  exceeded  that  of  the  ring  by  1/8  in. 
which  permitted  the  thermal  insulation  of  the  coil  with  a  uniform  layer 
of  a  fibrous  potassium  titanate  produced  by  Dupont  under  the  trade  name 
of  Typersul. 

The  coils  functioned  as  the  primary  of  a  transformer  and  the 
model  as  a  loaded  secondary.  This  heated  the  outer  perimeter  of  the 
rings  to  a  depth  of  0.030  in.  The  power  source  was  a  5  KW  radio  fre- 
quency generator  with  an  operating  frequenc/  of  0.375  MC. 

The  cooling  fixture  consisted  of  a  hollow  aluminum  hub,  inter- 
nally cooled  by  water,  to  which  the  interior  circumference  of  the  rings 
were  attached.  The  flow  of  water  through  the  hub  could  be  controlled. 
Preliminary  tests,  with  volumes  of  water  that  varied  between  0.2  and  2.0 
gal/min,  showed  that  the  amount  of  flow  had  little  influence  on  the 
temperature  gradient  that  could  be  maintained.  This  we  believe  was  due 
to  the  fact  that  the  hub  together  with  its  stand  were  made  of  aluminum, 
thus  providing  a  large  area  for  effective  heat  dissipation  to  the  air. 
Figure  30  is  a  close-up  view  of  the  model,  induction  heating  coils, 
insulator  and  cooling  system. 

d)  The  control  system 

Iron-constant an  thermocouples  were  fixed  by  crimping  to  the 
lower  face  of  the  rings.  The  cold  joint  was  kept  at  32   F.  All  thermo- 
couple wires  were  connected  through  a  radio  frequency  filter  and  a 
selection  switch  to  a  Leeds  and  Northroup  bridge  balance. 

In  every  ring  five  thermocouples  were  attached  along  a  radius. 
One  was  placed  at  the  outer  radius,  one  at  the  interior  radius,  and  the 
remaining  three  at  regular  intervals  between  the  first  two. 


6l4 


Errors  involved  in  the  determination  of  the  temperature  distri- 
bution were  less  than  1%.     The  steady  state  temperature  distribution  of 
Ring  No.  1  and  Ring  No.  2  were  compared  with  the  theoretical  distribution, 
xinder  the  assumption  that  the  coefficient  of  thennal  conductivity  is 
independent  of  the  temperature.  Good  correlations  were  obtained  as  is 
shown  in  Figures  31  and  32. . 

To  control  the  uniformity  of  the  heat  input  around  the  outer 
circumference  and  the  output  around  the  interior  circumference,  some 
rings  were  prepared  with  thermocouples  fixed  on  the  external  and  inter- 
nal perimeter  separated  ninety  degrees  of  arc  one  from  the  other. 
After  steady  state  heat  conduction  had  been  achieved  the  thermocouple 
readings  indicated  that  the  temperature  at  the  outer  circumference  was 
uniform  to  within  +3.$%   and  the  interior  circumference  to  within  +  k.O^- 

It.  2  The  Operation  of  the  Test 

For  each  ring,  the  radio  frequency  generator  was  set  at  2  KW  to 
begin  heating.  The  temperatures  were  read  off  the  thennocouples  and 
sufficient  time  was  given  for  the  gradient  to  stabilize .  Once  the 
temperatures  remained  stationary  over  a  period  of  fifteen  minutes, 
steady  state  conditions  were  assumed  to  have  been  obtained.  A  picture 
of  the  fringe  pattern  was  then  taken  and  the  temperatures  recorded. 

The  plate  voltage  of  the  generator  was  then  increased  by  0.5  KW. 
Again  time  was  given  for  the  thermal  circuit  to  reach  a  steady  state 
condition  and  another  picture  was  taken,  recording  the  corresponding 
temperatures . 
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k.3     Theoretical  Solutions  to  the 

Thermal  Stress  Problem 

The  interest  of  theoreticians  in  finding  a  solution  to  this 
thermal  botindaiy  value  problem,  considering  the  material  nonhoraogeneous 
due  to  the  temperature  and  therefore  radial  dependence  of  E,  v  and  a, 
has  been  considerable.  However,  experimental  data  to  judge  the  relative 
merit  of  the  necessary  simplifying  assumptions  has  been  lacking. 

Hilton  [38]  has  obtained  a  closed  form  solution  for  a  hollow 
cylinder  subjected  to  a  steady  state  temperature  difference  between  the 
interior  and  exterior  wall  under  the  assumption  that  E  and  a   are  arbi- 
trary functions  of  the  temperature  and  v  =  0.5-  He  also  presented 
a  series  solution  for  the  corresponding  plane  stress  problem. 

Chang  and  Chu  [39]   studied  independently  essentially  the  same 
problem  with  both  mechanical  as  well  as  thermal  loading.  They  did  not 
restrict  Poisson's  ratio  to  any  particular  value.  However,  the  tempera- 
ture distribution  they  assumed  across  the  thickness  of  the  cylinder  wall 
does  not  appear  to  be  physically  meaningful  if  steady  state  heat  trans- 
fer is  supposed  to  exist. 

Trostel  [I4O]  solved  the  problem  of  thermal  stresses  in  a  thick- 
walled  pipe  for  which  v  =  0.5j  the  coefficient  of  thermal  conductivity 
and  thermal  expansion  are  linear  fxinctions  of  the  temperature,  and 
Young's  modulus  is  a  quadratic  function  of  the  temperature. 

Nowacki  [I4I]  gives  further  references  in  the  Polish  and  Russian 
literature. 

Our  temperature  measurements,  Figures  31  and  3^,    agree  reason- 
ably well  with  the  solution  to  the  heat  conduction  equation  under  the 
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assumption  that  the  coefficient  of  thermal  conductivity  is  independent 
of  the  temperature  in  the  range  considered.  This  physical  fact  was 
assumed  correctly  by  Hilton  and  therefore  our  experimental  values  were 
compared  with  those  predicted  by  his  theory. 

This  comparison  hovrever  is  of  a  qualitative  nature,  since 
Hilton's  assumption  that  the  material  is  incompressible  does  not  pennit 
a  quantitative  comparison  betvreen  the  theoretical  values  of  the  stresses 
and  those  measured  in  a  real  material. 

All  investigators  agree  that  the  effect  of  a  temperature  depend- 
ent G  and  en   is  to  produce  different  values  of  the  stresses  across  the 
ring  than  those  the  elastic  homogeneous  theory  with  constant  G  and  a 
would  predict. 

Whether  the  stresses  increase  or  decrease  will  depend  on  the 
relative  variations  of  G  and  o'  with  temperature.  No  definite  statement 
can  be  made,  one  way  or  the  other,  unless  particular  numerical  values 
are  considered. 

For  both  rings  Hilton's  plain-strain  solution  was  computed, 
assuming  a  linear  variation  for  G  and  ex,   and  v  =  0.5.  Although  in  this 
case  the  plane-strain  solution  and  the  plane-stress  solution  will  give 
different  results,  numerical  computation  of  the  plane-stress  solution 
was  not  attempted  in  view  of  the  lengthy  calculations  involved  and  the 
poor  rate  of  convergence  that  Hilton  has  disclosed  in  his  paper. 

These  values  were  compared  with  the  elastic  solution  [Ii.2] 

assuming  G  =  G    equal  to  a  constant.  G    is  taken  as  the  shear 
°  avr  avr 

modulus  corresponding  to  the  average  temperature  in  the  cylinder. 
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The  coefficient  of  thermal  expansion  was  also  taken  as  ot  =  a 
using  for  the  average  value  the  same  definition  as  for  the  shear 
modulus.   V  was  taken  equal  to  0.5. 

For  Ring  No.  1  both  approaches  yielded  almost  identical  results 
•with  a  maximum  departure  of  h-3%-     The  results  corresponding  to 
Ring  No.  2   have  been  plotted  in  Figure  hO. 

The  value  of  the  yield  stress  as  a  function  of  temperature ,  for 
30k   stainless  steel,  is  as  follows: 

TABLE  3 

YIELD  STRESS  OF  30k   STAINLESS  STEEL  AS  A 
FUNCTION  OF  TEMPERATURE 

T  (°F)  a     (psi) 

7$  33,000 

500  31,500 

700  30,000 

900  2U,000 

1100  13,000 

1300  11,500 

1500  10,000 

The  values  in  this  table  must  be  kept  in  mind  -vdien  considering 
Figures  38,  39,   and  UO.  Wherever  the  thermal  stresses  exceed  the 
value  of  the  yield  stress  at  that  temperature,  the  numerical  value  of 
the  stress  has  no  significance  beyond  shoi-rf.ng  the  comparative  effect 
of  the  temperature  dependence  of  E,  v  and  a   [38]. 

We  must  also  realiae  that  the  values  reported  in  Table  3  have 
been  obtained  under  physical  conditions  xrtiich  are  different  from  those 
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prevailing  in  the  ring,  since  the  stresses  were  produced  by  mechanical 
loads.  Furthermore,  the  reported  yield  stresses  correspond  to  an 

arbitrary  definition  which  does  not  reflect  a  real  discontinuous  transi- 
tion in  the  behavior  of  the  material, 

U.i^  Experimental  Determination  of 
Thermal  Strains  and  Stresses 

The  moire  patterns  at  the  corresponding  temperatures  can  be 
seen  in  Figures  33   and  3U.  Degradation  of  the  image  due  to  the  process 
of  reproduction  is  evident. 

From  the  moire  fringe  patterns,  the  strains  shown  in  Figures  35 
and  36  were  computed,  using  the  continuous  light  intensity-displacement 
law  presented  in  [3]  and  [h] ■     The  details  of  the  automatic  processing 
of  the  information  were  presented  in  [6]. 

For  Ring  No.  1  the  strains  measured  experimentally.  Figure  3$, 
show  very  small  departures  from  the  corresponding  elastic  strains. 
This  is  interesting  in  view  of  the  fact  that  the  value  of  the  stresses 
would  indicate  that  plasticity  has  occurred  in  the  vicinity  of  the  inner 
and  outer  circumference.  Thus,  no  such  elastic  theory  is  rigorously 
applicable. 

For  Ring  No.  2  the  measured  strains,  Figure  36,  show  a  more 
definite  depsirture  from  those  predicted  by  the  elastic  theory  with 
constant  coefficients.  However,  the  elastic  theory,  vAiile  clearly  not 
applicable,  would  again  seem  to  yield  similar  values  for  the  strains. 

The  variation  of  the  mechanical  properties  of  the  stainless 
steel  AISI  type  30k   (l8Cr  -  8  Ni)  used  in  our  tests,  are  shown  in  Fig- 
ure 37.  These  values  were  obtained  from  references  [h3]   and  [hh] ' 
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Figures  38  and  39  show  the  comparison  of  theoretical  elastic 
stresses  with  constant  coefficients  with  stresses  determined  from  experi- 
mentally measured  strains  by  the  Dijhamel-Neumann  law  with  variables 
E,  V  and  O'. 

For  Ring  No.  1  the  differences  are  small.  This  coincides  with 
Hilton's  theory.  For  Ring  No.  2,   the  experimental  circumferential 
stresses  o     are  larger  than  those  of  the  homogeneous  elastic  theory  as 
predicted  by  Hilton's  theory.  The  radial  stresses  however  do  not  agree. 
Hilton's  solution  shows  an  increase  in  o     v^ile  our  experiment  indicates 
a  decrease  in  this  stress. 

In  interpreting  the  above  results^  the  reservations  made  in 
Section  5-3  concerning  the  meaning  of  the  stresses  when  the  yield  limit 
has  been  exceeded,  should  be  kept  in  mind.  In  Hilton's  numerical 
examples  the  stresses  are  also  of  such  magnitude  as  to  exceed  locally 
the  yield  stress  of  the  material. 

The  tests  reported  in  this  chapter  seem  to  indicate  that  for 
30ii  stainless  steel  up  to  1600  F  of  maximum  temperature  and  for  gra- 
dients of  680  F  or  less  the  value  of  the  thermal  strains  do  not  differ 
greatly  from  those  calculated  by  the  elastic  theory  using  E   ,  v 

and  a 

avr 

h.g  Residual  Strain  and  Stress  Determination 

Due  to  the  large  circumferential  stresses  present  at  both  the 
interior  and  exterior  perimeter  of  Ring  No.  2,  the  final  moire  pattern 
at  room  temperature  after  heating  was  analyssd. 
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The  resulting  strains  have  been  plotted  in  Figure  Ul.  The 
values  of  these  total  residual  circumferential  and  radial  strains  in 
the  vicinity  of  the  hole  and  the  outer  boundary  are   different  from 
zero  and  coincide  in  sign  with  the  strains  produced  by  the  thermal  load. 

The  strains  measured  in  this  manner,  due  to  the  geometrical 
character  of  the  moire  phenomenon,  are  the  final  resulting  strains  in 
the  body  that  was  initially  strain  free  before  the  loading.  As  such, 
they  are  the  sum  of  both  plastic  and  residual  elastic  strains  in  the 
body. 

In  order  to  separate  both  components  of  strain  the  retrieval 
of  the  elastic  strain  was  undertaken  by  means  of  a  destructive  test 
proposed  by  Sachs  [ii5]. 

This  test  consists  of  boring  out  the  central  portion  of  the  ring 
around  the  hole  and  measuring  the  variation  of  the  exterior  diaraeter 
for  a  given  bore  size. 

To  accomplish  this,  the  ring  was  placed  on  a  vertical  milling 
machine.  Figure  h^ ,    and  secured  at  the  outer  perimeter.  The  boring 
was  made  with  sharp  carbide  tool  bits  with  slow  machine  feeds  and  abun- 
dant cooling.  Two  jets  of  water  mixed  with  compressed  air  were  directed 
towards  the  interior  circumference  of  the  ring  for  heat  removal. 

After  each  cut  the  sample  was  allowed  to  cool  to  room  tempera- 
ture before  performing  a  measurement.  The  ambient  temperature  was 
recorded  to  withint  1/2  F  and  the  measured  values  corrected  accordingly. 

The  variations  in  the  external  diameter  were  measured  by  means 
of  a  Brown  and  Sharpe  transistorized  electronic  gage  capable  of  reading 
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0.000010  in.  The  instrument  was  connected  so  that  the  measurements 
would  not  be  affected  by  rigid  body  motions  of  the  test  sample  as 
a  whole. 

The  diametral  readings  were  taken  every  30  around  the  outer 
circumference  of  the  ring  and  averaged.  The  stresses  were  then  calcu- 
lated in  a  straight  forward  manner  [hS] .     The  plot  of  circumferential 
and  radial  stresses  is  shown  in  Figure  li3.  The  experimental  values  of 
the  circumferential  stress  satisfied  well  the  equilibrium  condition. 
From  the  stresses  the  residual  elastic  strains  were  computed.  Figure  Uli. 

Having  obtained  the  total  and  the  elastic  residual  strains,  the 
residual  plastic  strains  were  computed  according  to  the  expression 

T    E    P 
e  =  e  +  e 

The  residual  plastic  strains  are  shown  in  Figure  It5. 

The  analysis  of  Figures Ul,  k3,   hh   and  hS   show  that  the  total  and 
plastic  strains,  where  they  are  different  from  zero,  are  of  the  same 
sign  as  those  produced  by  the  thermal  gradient  applied.  The  elastic 
strains  and  stresses  are  of  opposite  sign  to  the  ones  produced  by  the 
thermal  load. 

The  ratio  of  plastic  circumferential  to  plastic  radial  strain 
at  the  edges  of  the  disk,  where  a  state  of  uniaxial  loading  can  be 
assumed  to  exist,  is  0.U50  inside  and  O.U?^  outside.  These  ratios 
differ  from  the  value  0.500  we  would  have  expected  to  obtain  under  the 
assumption  of  incompressible  plastic  flow. 
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The  reasons  for  this,  we  suspect,  are  the  accumulation  of 
experimental  errors  in  the  rather  indirect  process  by  vrfiich  the  plastic 
strains  are  found  and,  the  extrapolation  necessary  to  obtain  the  resid- 
ual, total  and  elastic  circumferential  strains  at  the  outer  edge  of  the 
ring  [1x6]. 


CHAPTER  V 
CONCLUSIONS  AND  RECOMMENDATIONS 

5-1  Conclusions 

In  this  research  the  analysis,  based  on  the  MTF  of  the  photo- 
optical  system  used  to  produce  and  record  the  moire  fringe  pattern, 
has  been  introduced  as  a  first  step  in  the  application  of  moire  to 
experimental  stress  analysis. 

It  has  resulted  in  a  better  understanding  of  the  optical  filter- 
ing capabilities  inherent  to  the  process,  and  it  provides  a  quantitative 
measure  of  the  limit  of  resolution  of  the  lens,  the  film  and  the  lens- 
film  combination.  Furthermore,  it  pennits  a  priori  information  on  the 
amount  of  contrast  that  can  be  expected  in  the  image  of  grids  or  fringes 
of  a  given  spatial  wavelength. 

The  strong  dependence  of  the  MTF  on  the  proper  focus,  the  paral- 
lelism of  object  and  image  plane,  the  rigidity  of  the  system  and  the 
aperture,  were  demonstrated  numerically. 

The  comparison  of  the  performances  of  the  ideal  and  real  lens- 
film  combinations  reveals  large  differences  that  amply  justify  the  extra 
effort  of  optical  testing  prior  to  a  moire  application. 

This  is  particularly  Important  in  problems  of  stress  concentra- 
tion where  the  displacement  gradient  around  notches  or  cracks  can  become 
large.  Under  these  conditions,  the  distance  between  moire  fringes  may 
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be  so  small  that  their  contrast  in  the  image  plane  can  be  drastically- 
reduced. 

Though  initial  mismatches  of  adequate  sign  could  correct  this 
problem,  it  is  clear  that  for  a  given  pitch  and  no  mismatch,  there 
exists  for  every  photo-optical  system  a  limit  to  the  displacement  that 
can  be  measured  dependent  on  the  MTF. 

In  this  work  the  first  steps  towards  effective  optical  filter- 
ing of  a  moire  signal  have  been  taken.  Again  the  MTF  appears  to  be  the 
logical  tool  for  the  evaluation  and  comparison  of  different  techniques. 

The  remarkable  properties  of  a  lens  in  a  coherent  light  field, 
as  an  analog  device  that  performs  a  Fourier  transform  and  displays  the 
Fourier  spectrum  of  a  signal,  has  been  brought  to  the  attention  of 
experimental  stress  analysts. 

This  idea  has  been  applied  with  success  to  the  filtering  of  the 
grid  lines  in  a  moire  fringe  pattern. 

An  incoherent  light  apodization  filter  has  been  implemented 
empirically  to  retrieve  periodic  line  patterns  immersed  in  strong 
random  noise.  The  potential  of  this  filtering  technique  was  thus 
brought  to  bear  on  a  problem  of  particular  interest  in  moire  work. 

Methods  of  numerical  in-phase  and  in-quadrature  filtering  were 
applied  to  overcome  certain  shortcomings  of  previously  used  electronic 
in-phase  filters.  With  them  the  phase  information  contained  in  the 
moire  light  intensity  trace  has  been  obtained  without  regard  to 
amplitude . 
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As  a  result,  the  minimum  phase  angle  that  can  be  meaningfully 

measured  has  been  reduced  from  5  in  previous  work  [6]  to  1  .  This 
implies  that  we  have  effectively  divided  the  grid  pitch  into  360  equal 
parts. 

The  accuracy  of  the  displacement  determined  at  this  level  of 
interpolation  for  a  carefully  tuned  filter,  has  been  found  to  be  no 
more  than  30  x  10   p  and  no  less  than  100  x  10   p. 

The  process  of  differentiation  by  which  strains  are  obtained 
from  the  experimental  curve  of  displacements,  was  critically  reviewed. 
The  solutions  presented  offer  a  degree  of  accuracy  and  reproducibility 
in  excess  of  previous  results  in  moire  work.  An  average  error  of 
50  X  10"  in. /in.  was  determined  for  the  strains  in  an  accurately  tuned 
filter. 

The  thermal  strain  problem  of  a  ring  under  axisymmetrical  steady 
state  temperature  gradient  was  solved  completely  by  the  moire  method. 
Interesting  results  on  the  similarity  in  the  magnitude  of  the  strains 
predicted  by  the  linear,  homogeneous  and  isotropic  theory  of  elasticity 
and  their  expeirimentally  determined  values,  have  been  reported. 

Since  the  stresses  in  the  ring  during  the  application  of  the 
thermal  load  exceeded  at  some  points  the  value  of  the  yield  stress  at 
the  corresponding  temperature,  the  total  residual  strains  were  obtained 
from  the  final  moire  pattern.  The  residual  elastic  strains  were 
obtained  by  the  boring-out  method  and  the  residual  plastic  strains  were 
then  obtained  by  superposition. 
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This  experiment  shows  that  when  elastic  and  plastic  components 
of   deformation  coexist  in  a  body,  moire',  due  to  its  geometrical  charac 
ter,  will  yield  the  algebraic  sum  of  the  strains.  The  separation  into 
their  respective  components  can  be  accomplished  only  after  one  of  the 
two  is  measured  by  means  other  than  the  nondestructive  analysis  of  the 
final  moire  pattern. 

The  techniques  of  model  grid  engravure  able  to  withstand  high 
temperatures,  have  been  presented  in  detail.  The  photoengraving  tech- 
niques, the  composition  of  etchant  baths,  the  engravure  times  and  the 
current  densities  reported  have  produced  consistently  good  results. 
We  hope  they  will  provide  useful  guidelines  for  workers  in  this  field. 

5.2  Recommendations 

The  ability  of  moire  to  measure  displacements  is  well  known  [h7, 
hQ,  h9],   however,  its  application  to  strain  analysis  implies  a  process 
of  differentiation  which  is  prone  to  error. 

The  methods  of  hand  and  numerical  derivation  are  costly  in  real 
and  computer  time,  respectively,  therefore,  efforts  should  be  made  to 
implement  optically  this  procedure. 

One  such  technique  was  proposed  by  Dantu  [50]  and  another  by 
Parks  and  Durelli  [5l].  Both,  however,  involve  a  rigid  body  motion  of 
the  moire'  pattern  or  the  deformed  grid  with  respect  to  itself  to  pro- 
duce fringes  that  are  the  loci  of  points  of  equal  strain.   This  shift 
results  in  a  slight  indetermination  as  to  what  point  in  the  body  is  sub- 
jected to  a  given  strain. 
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An  optical  method  of  derivation  that  does  not  require  this  dis- 
placement would  be  ideal.  In  this  connection,  the  optical  properties  of 
coherent  light  systems  may  yield  the  basis  for  an  alternate  solution. 

The  minimum  photographic,  optical,  electronic  and  digital 
instrumentation  necessary  to  produce,  record  and  analyse  moire'  patterns 
along  the  lines  proposed  in  this  research  should  be  determined.  The 
apparatus  used  here  is  considered  to  exceed  in  many  cases  what  is  neces- 
sary to  obtain  results  of  engineering  accuracy. 

In  the  numerical  retrieval  of  phases  from  the  moire'  light  inten- 
sity trace  it  is  conceivable  that  the  decomposition  of  the  signal  in 
other  than  a  trigonometric  series  could  yield  faster  convergence  and 
more  accurate  results.  This  matter  would  merit  further  analysis. 

Electronic  in-phase  and  in-quadrature  filters  should  be  applied 
for  the  processing  of  moire'  data  in  order  to  increase  the  versatility 
and  speed  of  the  analysis. 

If  the  statistical  properties  of  the  noise  in  the  moire'  light 
intensity  trace  were  known,  the  requirement  that  the  noise  and  signal 
spectrum  should  be  disjoint  could  be  relaxed.  This  would  permit  the 
implementation  of  a  filtering  technique  more  powerful  than  the  ones 
presented  here.  However,  the  task  of  characterizing  the  noise  content 
of  moire  fringe  patterns  obtained  under  diverse  conditions  of  practical 
interest  has  not  been  done  to  date. 

Based  on  the  complete  strain  determination  made  for  the  stain- 
less steel  ring,  a  thermoplastic  analysis  in  terms  of  stresses  can  be 
attempted.  This  will  require  the  knowledge  of  the  stress-strain  rela- 
tion in  the  elastic  and  plastic  region  of  behavior  of  the  metal. 
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Some  of  these  ideas  are  being  pursued  in  the  continuing  research 
effort  of  which  this  dissertation  is  a  part;  others  remain  as  avenues 
possibly  wirth  exploring  by  anyone  attracted  to  this  matter  of  large 
field  strain  analysis  by  a  method  of  such  direct  visual  appeal  as  moire. 


APPENDIX  A 

ON  LINEAR  FILTER  THEORY 

A.l     Linear  Systems 

For  the  purpose  of  this  appendix  a  linear  system  is  a  linear 
operator  L,    such  that  for  all  functions  cp(x)   and  \(r(x)   in  the   space 
domain  and  all  scalars  a,  b,   the  following  relation  holds 

L[a9(x)   +  bKx)]   =  aL[cp(x)]   =  +  bL[Kt)]  (A.l) 

and  that,   if 

L[Kx)]    =  cp(x)  (A  ^2) 

and  x^  is  a  real  constant,   then 

UHx  -  x^)]    =cp(x  -  x^)  (A  .3) 

that  is,  L  is  invariant  under  a  space  transformation. 
Then,  if 

L[6(x)]  =  h(x)  (A.i,) 

L[6(x-x^)]  =h(x-x^)  (A.5) 

it  follows  that 


cp(x)  =  r   Ke)h(x  -  e)  dc 


_00 


(A. 6) 


If  we  call  cp(x)  the  output  of  the  linear  system  L  and  Hx) 
the  input,  (A  .6)  indicates  that  the  output  of  L  is  given  in  terras 
of  the  input  and  a  unique  function  h(x)  called  the  weight  function 
or  impulse  response  of  L. 
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The  Fourier  transform  of  h(x)  is  called  the  transfer  function 

of  L  and  is  given  by- 
es 
H(f)  =  £„  h(x)e-2"i^^  dx  (A. 7) 

If  we  recognize  (A. 6)  as  the  convolution  product  we  can  write 

cp(x)  =  Kx)  ^  h(x)  (A. 8) 

and 

f(f)  =  Y(f)  •  H(f)  (A.9) 

i*ere  $(f),  "^(f)  and  H(f)  are  the  Fourier  transforms  of  9(x), 
Kx)  and  h(x)  in  that  order. 

From  (A.9)  we  obtain  immediately  the  transfer  function  that 
will  produce  a  desired  output  whose  transform  is  $(f)  from  an  input 
Hf)   is 

H(f)  =  M-  (A.IO) 


The  expression  for  H(f)  is  in  general  complex  and  of  the  form 

H(f)  =  A(f)e^^^^^  (A. 11) 

We  will  restrict  our  attention  to  linear  transformations 

involving  no  phase  shift,  thus 

e(f)  =  0  (A  .12) 

We  define  a  low-pass  filter  as  a  filter  for  vrfiich  A(f)   is 

small  in  some  sense  for   If     >  f   ,   where  f     is  called  the  cut-off 

c'  c 

frequency.     A  band-pass  filter  is  a  filter  where  A(f)  is   small  outside 
an  interval    |f  +  -s-| ,  Figure  l8. 
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The  ideal  band-pass  filter  would  have  the  follovd.ng  transfer 
function.  Figure  3-la 


H(f)  =  A(f) 


'i  If„-^l5kl<k,^^l 


c  U^.^|>|f|>U^.f| 


(A.13) 


A. 2  Definition  of  a  Digital  Filter 

If  the  function  Kx)  defines  a  generalized  function,  is  band 
limited  and  we  possess  its  sampled  values  obtained  at  a  sampling  fre- 
quency f  which  is  at  least  twice  as  large  as  the  largest  frequency  f 
present  in  the  function,  Shannon's  sampling  theorem  guarsintees  that  '^{x) 
can  be  uniquely  determined  from  its  sampled  values 

*n  ==  ^^2^)  =  ^^f) 
a  s 

by  the  Infinite  series 

Kx)  =  "  „  a  *  (A.lii) 

n=  _«>   n  n 

If  H(f)  is  such  that  ^(f)    can  be  written  as  a  trigonometric 
series,  then 

2nTTi^ 
H(f)  =5]    a„  e  (A. 15) 


n=  -00 


n 


The  transform  of  (A  .15)  will  be 


h(x)  =  E    a  6(x  +  :^)  (A  .16) 

n=  -00   n       X 
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If  we  define  the  convolution  h(x)  ■«•  i|f  ,  where  i  is  the  space 

sampled  version  of  K^),  then  the  sampled  version  of  the  desired  output 

9  =  9(J^)  will  be  given  by 
s 

cp(f )  =  "    a  K^)  (A.17) 

J-g    n=  -«>      -^s 

For  computational  purposes  the  infinite  series  (A  .15) ,  (A. 16) 
and  (A.17)  have  to  be  replaced  by  finite  summations,  then  (A  .15) 
becomes  „ 

H„(f)  =  E   a^  e     ^  (A. 18) 


n=-N  " 


and  (A.17) 


N 

^  =  S   a  1(1  (A  .19) 

n   n=  -N  "^  "i+J^ 

where  cp  is  a  good  approximation  of  9(-»-)  through  (A. 19). 

s 
Equation  (A. 19)  is  the  fundamental  formula  of  digital  filter- 
ing [33]  and  is  sometimes  referred  to  as  a  finite  moving  average 
process  [32]. 


APPENDIX  B 


ON  THE  PREPARATION  OF  METAL  MODELS 
FOR  HIGH  TEMPERATURE  MOIRJ!  TESTS 


B . 1  Introduc  tl on 

Whatever  means  we  have  for  analyzing  the  moire  data,  we  are 
interested  in  obtaining  the  fringe  patterns  as  free  as  possible  of  noise. 
This  noise  obscures  the  meaningful  signal  and  makes  the  processing  of 
th^  information  more  difficult  and  less  accurate. 

If,  as  is  customary,  a  photograph  is  used  as  a  means  of  memo- 
rizing the  moire  pattern,  this  noise  will  appear  as  spots,  blemishes 
or  variations  of  background  intensity  in  the  negative. 

High  temperature  tests  are  especially  prone  to  this  inconveniency 
due  to  the  fact  that  organic  and  inorganic  substances  used  to  print  lines 
on  the  model  change  color  with  temperature. 

The  non-uniformity  of  the  temperature  fields  investigated  will 
cause  a  variation  of  the  model's  color  that  will  show  as  low  spatial 
frequency  variations  of  the  background  intensity  on  the  photograph.  Dust, 
scratches  or  pits  on  the  model  surface  will  appear  as  high  frequency 
disturbances  on  the  recorded  pattern  of  fringes. 

To  obtain  the  best  moire,  over  a  wide  range  of  temperatures, 
a  heat  resistant  grid  of  excellent  quality  must  be  engraved  on  the 
model.  This  involves  the  careful  completion  of  the  following  three  steps: 
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Surface  finishing  of  the  model 
Photo-engravure  of  the  grid 
Chemical  etching  of  the  grid 

B.2  Surface  Finishing  of  the  Model 

In  general  it  can  be  said  that  the  surface  of  the  model  that  is 
to  be  engraved  with  lines  should  be  polished  until  free  from  pits  and 
scratches.  The  amount  of  remaining  unevenness  that  can  be  tolerated 
will  depend  on  how  this  shows  on  the  picture  of  the  moire  pattern.  Also, 
if  the  surface  roughness  hinders  the  printing  of  the  lines,  or  if  it 
permits  under  etching  of  the  photo-resist,  finer  polishing  will  be 
called  for. 

A  decrease  in  the  grid  pitch  will  necessitate  an  increase  in 
the  surface  polish.  We  can  appreciate,  here  again,  the  advantage  of 
being  able  to  work  with  fewer  lines  per  inch  for  a  given  sensitivity. 

<    Softer  metals,  like  copper  and  aluminum,  can  be  finished  by- 
hand  on  a  lapping  table  if  the  model  is  not  too  large.  The  usual  pro- 
cedure is  to  polish  them  with  progressively  finer  emery  paper  until  we 
obtain  a  suitable  finish.  For  300  lines  per  inch  this  was  found  to  be 
grit  No.  hOO  A. 

Harder  metals  and  samples  of  larger  size  cannot  be  easily 
polished  by  hand.  The  stainless  steel  models  prepared  for  the  tests 
in  the  main  text  were  ground  and  polished  mechanically  to  a  15  micro 
inch  finish  which  proved  sufficient  for  engraving  300  lines  per  inch. 

It  was  found  necessary  to  nickel  plate  electrically  or  chem- 
ically all  mild  steel  test  pieces  for  thermal  testing.  This  prevented 
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the  oxidation  of  their  surface  that  reduced  the  contrast  between  the 
grid  lines.  The  thin  plated  layer  was  then  polished  and  engraved. 

For  more  critical  engraving  work,  such  as  750  lines  per  inch 
in  metal  samples,  it  was  necessary  to  increase  the  mechanical  grinding 
and  polishing  to  obtain  a  6  micro  inch  finish  for  good  engravings. 
A  typical  example  of  insufficient  surface  polishing  in  an  aluminum 
sample  engraved  with  a  cross  grid  of  750  lines  per  inch  can  be  seen  in 
Figure  ii6. 

A  very  satisfactory  surface  finish  can  be  obtained  by  electro- 
polishing  the  surface  to  be  engraved.  Almost  without  exception  etch- 
ing baths  will  serve  as  good  polishers  if  conveniently  diluted,  or  if 
the  current  densities  used  for  etching  are  reduced  to  an  adequate  level 
to  permit  polishing. 

Such  is  the  case  of  the  sodium  hydroxide  bath  for  aluminum  and 
the  orthophosphoric  bath  for  stainless  steel  desciribed  below. 
Aluminum  polishing: 

2  oz  sodium  hydroxide  pellets 
1  gal  water 
30  amp  per  sq  ft  current  density 
Stainless  steel  polishing: 

1  vol  orthophosphoric  acid  {8S%   cone) 
6  vol  water 
72  amp  per  sq  ft  current  density 
If  the  polishing  times  are  found  to  be  too  long  for  some  appli- 
cations, the  bath  can  be  activated  either  electrically  or  by  increasing 
the  acid  or  alkali  concentration. 
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The  work  to  be  polished  should  present  a  uniform  surface  finish- 
ing, with  no  deep  scratches  or  pits,  before  immersion  in  the  bath,  since 
the  action  of  the  electropolishing  will  not  eliminate,  but  rather  accen- 
tuate, any  large  crevice  in  the  surface. 

By  careful  electropolishing  surfaces  can  be  finished  to  receive 
photo-engraved  or  chemically  milled  750  to  1000  lines  per  inch  by  the 
processes  that  will  now  be  described. 

B.3  Photo-Engravure  of  the  Grid 

Photo-resists  have  two  properties  that  make  them  useful  for  moire 
work.  One  is  their  ability  to  reproduce  photographically  by  contact  a 
grid  printed  on  a  negative  which  we  call  film  master.  Another  is  their 
resistance  to  acid  and  alkali  that  permits  us  to  etch  chemically  the 
areas  of  the  metal  surface  between  the  lines  protected  by  resist. 

TVo  commercially  available  photo-resists  have  been  employed  exten- 
sively t)y  us.   These  are  the  Kodak  Photo  Resist  manufactured  by  the 
Eastman  Kodak  Company  of  Rochester,  New  York,  and  the  Dynachem  Photo 
Resist,  distributed  by  Lithoplate,  Inc.,  of  Covina,  California,  for  the 
Dynachem  Corporation. 

Both  require  handling  in  subdued  tungsten  illumination.  A  1^ 
bulb  not  less  than  I   ft  away  proved  to  be  satisfactory. 

A  listing  of  operations  and  times  that  consistently  yielded 
very  good  photo-engravures  follows. 

Photo-engravure  with  Kodak  Photo  Resist  (KPR): 

a)  Spray  or  whirl  on  a  fine  even  coat  of  KPR  and  let  dry. 

b)  Expose  master  grid  and  model  in  vacuum  frame  to  20  amp 
carbon  arc  lamp  at  1  ft  distance  for  U  min. 
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c)  Develop  with  KPR  developer  for  at  least  h   mixi.  Longer 
developing  times  have  no  effect  on  the  image  quality. 

d)  Flow  with  KPR  dye  for  1  min. 

e)  Wash  with  water. 

f)  Blast  with  air  to  dry. 

g)  Heat  in  oven  for  30  min  at  li70°F.  Then  allow  to  cool 
to  room  temperature  before  further  handling. 

Photo-engravure  with  Dynachem  Corporation  Resist  (DCR): 

1   a)  Spray  the  surface  placed  vertically  with  a  slow  horizontal 
sweeping  motion,  advancing  from  bottom  to  top. 

b)  Lay  sample  flat  for  2  min  to  allow  the  dry-wet  coating  to 
settle  to  a  uniform  thickness. 

c)  Place  specimen  in  an  oven  or  on  an  open  electric  hot  plate 
for  8  to  10  rain  at  lIiO°F.  Then  allow  to  cool. 

d)  Expose  in  vacuum  frame  with  master  for  7  min. 

e)  Develop  in  DCR  developer  for  30  to  60  sec.  The  develop- 
ing time  is  critical. 

f)  Spray  on  DCR  dye.  Leave  dye  on  surface  for  5  sec. 

g)  Wash  in  water. 

h)  Blast  with  air  to  dry. 

i)  Place  in  an  oven  or  on  a  hot  plate  and  heat  for  30  min 
at  2ii0  F.  Then  allow  to  cool  to  room  temperature  before 
further  handling. 

For  both  the  KPR  and  DCR  processes  the  optimum  exposure  and 
developing  time  will  depend  on  the  thickness  of  the  coating  applied. 
Since  the  thickness  of  this  thin  film  is  difficult  to  monitor  by  ele- 
mentary means,  some  experimentation  will  always  be  necessary  to  obtain 
the  best  results. 
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Figure  U?  is  a  UOX  magnification  of  a  dyed  photo-engraved 
crossed  grid  of  300  lines  per  inch.  The  rounding  of  the  comers  is 
due  to  light  diffusion  in  the  photo-sensitive  coating  thickness. 

Some  manufacturers  of  photo-resists  recommend  that  the  resist 
should  not  be  dyed  if  maximum  resistance  to  the  etchant  is  desired. 
This  precaution  was  taken  in  the  engravure  of  the  stainless  steel  rings 
of  the  main  text.  The  evaluation  of  the  quality  of  the  photo-engraved 
grid  was  performed  by  inspection  of  the  moire  pattern  produced  by  the 
engravure  with  the  master.  The  usual  inspection  by  microscope  proved 
to  be  very  difficult  since  the  resist  is  colorless  and  transparent. 

In  an  earlier  paper  [^^] ,   vAiere  a  thermal  problem  of  the  same 
geometry  as  the  one  presented  in  this  dissertation  was  studied,  1000 
lines  per  inch  were  photo-engraved  with  dyed  KPR  as  explained  above. 
This  printing  process  produced  moire  patterns  up  to  UOO  F.  However, 
the  changes  in  color  of  the  organic  coating  did  produce  some  undesirable 
background  intensity  variations. 

Since  most  photo-resists  clear  at  about  600  F.  and  turn  to  white 
ash  at  7^0  F,  500  F  seems  to  be  the  limit  of  applicability  for  this' 
method  of  engraving  thermal  test  pieces. 

B.U  Chemical  Etching  of  the  Grid 

Having  completed  all  the  steps  described  in  Section  B.3,   the 
sample  can  be  immersed  in  an  acid  or  alkali  etch  bath. 

The  following  solutions  have  been  found  to  be  effective  etchants 
[53]: 
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For  aluminum,  a  }>%   concentration  caustic  soda  solution,  made 
by  dissolving  Ii  oz  sodium  hydroxide  pellets  in  1  gal  water.  For  a  more 
active  bath  the  concentration  can  be  increased  to  S%   by  dissolving 
12  oz  sodium  hydroxide  per  1  gal  water. 

For  copper,  an  8^  concentration  hydrochloric  acid,  prepared  by 
adding  10  oz  38^  hydrochloric  acid  to  1  gal  water. 

For  nickel  and  copper,  1  gal  photoengravers  iron  chloride 
(iiO  Be)  added  to  6  oz  38^  hydrochloric  acid.  Keep  the  temperature  of 
the  bath  below  IIO't  by  suitable  cooling. 

For  mild  steel  [5h-],  a  S%   concentration  nitric  acid  solution. 

For  stainless  steel,  a  solution  of  '}>'})%   orthophosphoric  acid, 
Sl%   of  specific  gravity  1.26,  and  10^  water.  A  current  density  of 
72  amp  per  sq  ft  is  applied  to  the  sample  made  anodic. 

The  time  the  model  can  be  left  in  the  bath  will  depend  on 

a)  the  depth  of  etching  required 

b)  the  resistance  of  the  photo  resist 

c )  the  amount  of  under-etching  that  develops . 

a)  The  bath,  according  to  its  concentration  and  the  current 
density  applied  to  the  sample,  removes  material  at  a  constant  rate, 
e.g.,  a  bath  of  1  volume  of  38^  hydrochloric  acid  and  U  volumes  of  water, 
at  \x   amps  per  sq  in.  and  2  volts,  removes  material  from  a  stainless  steel 
plate  at  0.001  inch  per  minute.  By  periodic  inspection  of  the  etched 
surface  under  the  microscope,  the  frequency  of  the  inspections  will 
depend  on  the  activity  of  the  bath,  we  are  able  to  control  the  process 
and  remove  the  sample  lAien  the  lines  are  etched  to  sufficient  depth. 
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If  this  time  is  recorded,  similar  samples  will  take  approximataLy 
equal  times  to  be  etched  the  same  amount  as  long  as  the  bath  does  not 
become  exhausted.  Aging  of  the  bath  vrill  result  in  longer  etching 
times  if  the  current  densities  are  kept  constant. 

b)  What  was  said  above  implied  a  resist  vrfiich  would  withstand 
the  action  of  the  bath  as  long  as  we  desired.  This  will  not  be  true  if 
the  cleanliness  prior  to  application  of  the  photo-resist  were  not  perfect, 
the  ^thickness  of  the  protective  layer  were   not  enough,  the  exposure 
incorrect,  the  baking  insufficient  or  abrasion  had  destroyed  the  contin- 
uity of  the  engraved  grid. 

It  is  important  to  note  that  imderexposure  of  the  photo-resist 
causes  the  developer  to  wash  away  protective  coating,  while  overexposure 
will  cause  stripping  of  the  coating  in  the  wash  after  developing. 

All  the  above  mentioned  causes  will  diminish  the  resistance  of 
the  photo-resist  and  may  cause  it  to  strip  off  the  surface  it  is  pro- 
tecting before  the  necessary  depth  of  etching  has  been  achieved. 

When  grids  are  being  etched  the  stripping  of  the  resist  in  some 
localized  areas  will  force  us  to  interrupt  the  etching  entirely.  Fail- 
ure to  do  this  will  cause  indiscriminate  etching  of  the  unprotected 
metal,  thus  erasing  whatever  line  pattern  had  been  etched  already. 
Longer  periods  of  etching  will  eventually  cause  pitting  in  the  unpro- 
tected areas  and  more  stripping  due  to  underetching  as  id.ll  be  explained. 

c)  We  have'  called  underetching,  for  lack  of  a  better  name,  the 
attack  on  the  sides  of  the  partially  etched  trough  by  the  etchant,  that 
causes  metal  adhered  to  the  resist  on  the  model  surface  to  be  removed 
from  below.  Figure  1^8  a  and  Ii8b. 
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The  photo-resist  can  then  strip  off  by  the  removal  of  metal 
to  which  It  was  adhered. 

We  have  found  in  repeated  testing  that  the  raaximum  depth  of 
etching  without  significant  underetching  for  Ronchi  type  rulings  varies 
between  0.5  p  for  KPR  to  0.75  p  for  DOR. 

Failure  to  extract  the  sample  from  the  etch  bath  when  under- 
etching  has  begun  will  produce  an  increase  of  the  width  of  the  etched 
lines  making  the  grid  lose  its  Ronchi  characteristic,  Figure  USc. 
Further  underetching  will  produce  rounding  off  of  the  line  edges  by 
gradual  stripping  of  the  underetched  resist.  This  will  diminish  the 
contrast  between  lines  of  the  engraved  grid  and  upon  superposition  with 
the  master,  a  less  well-defined  moire  will  result. 

As  we  have  seen,  the  depth  of  etching  depends  on  various  fac- 
tors all  of  which  must  be  carefully  controlled.  The  main  difficulty  in 
the  production  of  a  good,  deeply  etched  model  is  the  stripping  of  the 
resist  in  the  etch  bath  before  sufficient  etching  has  occurred. 

As  we  said  before,  even  if  the  stripping  is  local,  the  uniform- 
ity of  the  grid  will  be  affected  if  the  model  is  not  removed  immediately 
from  the  solution. 

If  etching  has  not  begun  before  the  resist  strips  off,  we  will 
be  a1)le  to  remove  the  remainder  completely,  and,  after  cleaning  thoroughly 
the  metal  surface,  proceed  to  repeat  the  photo-engravure  correcting 
■vriiatever  caused  the  resist  to  fail. 

However,  if  the  etching  had  begun  before  the  masking  failed, 
but  had  not  progressed  to  enough  depth  for  our  purpose,  we  are  left 
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with  a  shallow  etched  surface.  Since  this  etching  has  destroyed  our 
surface  finish,  we  are  obliged  to  return  our  sample  to  the  very  first 
stage  of  its  preparation.  This  not  uncommon  occurrence  is  dishearten- 
ing and  time  consuming. 

There  is  a  method  used  by  engravers  that  increases  the  resist- 
ance of  the  resist  during  the  etching  process  that  corrects  small  defi- 
ciencies that  might  develop  in  the  continuity  of  the  coating.  It 
involves  the  use  of  a  substance  known  in  the  trade  as  Dragon's  Blood. 
This  is  a  very  fine  organic  red  powder  that  can  be  applied  on  the  wet 
surface  that  is  being  etched  with  a  soft  camel's  hair  brush.  The  motion 
of  the  brush  should  be  perpendicular  to  the  direction  of  the  grid  lines 
allowing  the  powder  to  dissolve  and  be  swept  against  the  sides  of  the 
lines  and  into  any  small  pit  in  the  resist. 

If  the  model's  surface  is  then  dried  with  hot  air,  the  dis- 
solved powder  will  adhere  to  the  resist  and  build  up  on  the  etched  cor- 
ners, as  can  be  seen  in  Figure  USd.  This  added  protection  will  per*mit 
deeper  etching. 

Though  this  treatment  has  been  proven  effective  for  etching' 
coarse  grids  of  50  to  150  lines  per  inch,  its  use  has  not  yielded  sig- 
nificant improvement  for  denser  grids.  It  seems  obvious  that  when  the 
interline  spacing  is  reduced  it  becomes  increasingly  difficult  for  the 
dissolved  powder  to  get  between  the  lines. 

We  would  wish  to  be  able  to  take  a  shallow  etched  surface  and 
photo-engrave  a  new  resist  grid,  identical  to  the  one  that  stripped 
off,  and  continue  etching.  This  implies  the  necessity  of  placing  our 
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master  exactly  in  the  same  position  with  respect  to  the  model's  surface 
as  it  was  when  we  exposed  the  now  destroyed  grid. 

Because  the  master  is  made  of  straight  parallel  lines  of  constant 
pitch,  free  of  any  significant  periodic  error,  the  above  correspondence 
requires  only  the  exact  superposition  of  one  line  of  the  master  with  one 
line  of  the  semi-etched  model.  It  is  conceivable  that  such  superposition 
could  be  accomplished  with  the  aid  of  a  microscope  and  a  great  deal  of 
patience. 

A  macroscopic  method  of  superposition  has  been  employed  with 
success  for  this  same  purpose.  It  consists  of  superimposing  the  par- 
tially etched  model,  sensitized  with  unexposed  photo-resist,  with  a 
photographic  glass  copy  of  the  master  grid.  By  eliminating  all  moire 
fringes  due  to  rotation,  absolute  parallelism  of  both  sets  of  grids 
can  be  insured.  By  displacing  without  rotation  the  master  with  respect 
to  the  model,  a  change  in  the  light  intensity  reflected  from  the  model 
is  produced.  The  maximum  reflected  intensity  over  the  entire  model 
surface  will  indicate  the  position  in  which  the  correspondence  between 
the  partially  etched  lines  of  the  model  and  the  dark  lines  of  the 
master  is  perfect.  Upon  obtaining  this  correspondence  the  model  can  be 
exposed  in  intimate  contact  with  the  master  in  a  vacuiim  device.  The 
only  requirement  for  this  method  to  be  applicable  is  that  the  partial 
etching  be  sufficient  to  produce  at  least  a  faint  moire  pattern. 

When  only  shallow  etching  is  possible  in  a  metal  due  to  the  grid 
pitch  or  other  circumstances,  an  alternative  method  to  the  standard  pro- 
cedure described  above  can  be  used.  This  consists  of  selecting  the 
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reducing  or  oxidizing  atmosphere  and  temperature  that  will  cause  the 
metal  surface  to  change  color.  This  transformation  can  be  made  to  take 
place  only  to  a  very  shallow  depth  from  the  sample  surface  and  not 
destroy  the  surface  finish. 

Based  on  this  principle  polished  stainless  steel  samples  were 
introduced  in  a  hydrogen  furnace.  The  machine  used  was  a  Sargeant  and 
Wilbur  oven  manufactured  by  Sargeant  and  Wilbur,  Inc.,  of  Pantucket, 
Rhode  Island.  The  furnace  was  operated  at  1000  C  with  a  15  cu  ft  per 
rain  flow  of  wet  hydrogen.  The  samples  were  left  in  the  oven  for  30  min. 

This  treatment  left  the  surface  of  the  stainless  steel  colored 
dark  green.  The  depth  of  the  color  layer  depends  on  the  extent  of  the 
chrome  reduction  and  can  be  regulated  by  the  time  the  sample  is  left 
in  the  oven. 

On  removal  from  the  oven,  and  after  cooling,  the  samples  were 
photo-engraved  and  placed  in  the  etching  bath,  A  very  shallow  etch 
cut  through  the  thin  green  layer  exposing  the  bright  polished  metal 
below. 

The  final  step  in  the  preparation  of  samples  for  thermal  tests, 
after  the  etching  of  the  model  has  been  completed,  is  the  removal  of 
all  traces  of  resist  from  the  surface  prior  to  heating.  Failure  to  do 
this  will  result  in  the  organic  coating  first  changing  color  and  then 
burning  off  vrtiile  the  test  is  in  progress.  To  accomplish  this  cleaning 
KPR  "Strip-Ur"  SK  h   and  SK  5  distributed  by  Etchomatic,  Inc.,  for 
Spar  Kleen  Products,  Inc.,  was  found  to  be  effective. 
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Figure  h9   is  a  lOX  magnification  of  an  etched  cross  grid  of 
300  lines  per  inch  in  stainless  steel  with  all  traces  of  resist  removed. 
Some  underetching  has  caused  the  lines  to  narrow  between  crossings. 
The  rounded  appearance  of  the  etched  holes  is  due  to  the  diffraction  of 
light  during  photographic  printing  as  was  explained  above. 


APPENDIX  C 
PROGRAMMING  COEES 
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C 

c 
c 


GOODMAN  S  COMBINATION  OF  TUKEY  IN  PHASE  AND  IN 
QUADRATURE  FILTtRS 


C 
C 

c 
c 
c 
c 


c 

c 
c 
c 
c 
c 
c 
c 
c 
c 


BX(500),BY{500) ,DATA(4C00) ,TAK(40G0) 
FMT  (12),  TITLE  ( 7 )  .RAN { 40C0 )  ,CUM ( 100 ) 


DIMENSION 

DIMENSION 

PI=3.1416 

P8=8./PI 

P4  =  '!../PI 

RRD=1HT 


FMT=FORMAT  OF  DATA 
TITLE=riTLE  OF  DATA 
NDATA=NUMBER  UF  DATA 
DATA=DATA  POINT 

READ  INPUT  TAPE  5 , 25, ( FMT ( I)  ,  1  =  1  ,  12 ) 

25  FORMAT  ( 12A6) 

READ  INPUT  TAPE  5, 26, ( T 1 TLE ( I  )  , 1= 1 , 7  ) 

26  FORPAT  (7A6) 

READ  INPUT  TAPE  5,3,NDATA 
3  FORMATi 14) 

READ  INPUT  TAPE  5 , FMT , ( DATA ( I) , 1= I , NDATA ) 

RID=LETTER  T  IN  COLUMN  1  TO  IDENTIFY  TUKEY  S  PARAMETER 

CARD 

TERMS=  NUMBER  OF  TERMS  IN  TRIGONOMETRIC  EXPANSION  OF 

TRANSFER  FUNCTION 

CW=CENTER  WAVELENGTH  IN  DATA  POINTS 

NST=DEFINES  STARTING  DATA  POINT  NUMBER  FOR  ABOVE 

PARAMETERS 

NFN=DEFINES  END  DATA  POINT  NUMBER  FOR  ABOVE  PARAMETERS 

7  READ  INPUT  TAPE  5 , 1 , RI D , TERMS  ,0W ,NST  ,  NFN 
1  FORMATC A1,F9.0,F10.0,2I10) 
IF(RID-RRD)  13,14,13 

13  WRITE  OUTPUT  TAPE  6,15 

15  FORMATC lOX,  15H  LRROR  IN  NDATA) 
CALL  EXIT 

14  IF(NST+NFN)  11,12,11 
12  NST=1 

NFN=NDATA 
11  IF(TFRMS)  9,9,10 
10  M=TERMS-1. 

WRITE  OUTPUT  TAPE  6 , 27 , (  1 1 TLE { I ) , I  =  1  ,  7 ) 

27  FORMAT  ( IHl 7A6 , 6X 12HTUKt Y  FILTER//) 
WRITE  OUTPUT  TAPE  6, 26 , T ERM S , OW , NST, NFN 

28  FORMAT  ( 9H  TERMS  -  F10.6,6H  OW  =  F10.6,7H  NST  =  16, 
$7H  NFN  =  16///) 

WRITE  OUTPUT  TAPE  6,29 

29  FORMAT  ( X6HNUMBER6X lOHI NPUT  DATA13X7HVALUE  X15X 
i7HVALUE  Y14X 

$  9HRESULTANT15X5HTHETA//) 
W0=2.0*PI/0W 
PM=PI/TERMS 
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T8=8. /TERMS 

PM2=PM*2. 

TA=4. /TERMS 

CAST=-{0.5A-0.46)»T4 

TW=TERMS»WO 

CASTY=CAST*SIiNF(TW) 

CASTX=CAST«COSF(TW) 

AY  =  0. 

AX  =  0. 

CO    2    J=1,M 

R=J 

DD=R»PM 

D=CCSF(DD) 

C=0. 54  +  0. 46*0 

E=R*wa 

EY=SINF(E) 
EX=COSF(E) 
CO=C«D 

AY=AY+C»EY»EY*SINF(PM2»R)/R 
AX=AX+CC»£X»EX»SINF{DD)/R 
BY{ J)=CD»EY 
BX{ J)=CO»EX 
AY=AY»P4 
AY=1./AY 
AX=(P8»AX)+T4 
AX=1./AX 

FORMAT!  10X,E20.ci) 
TOTAL=0. 
RDATA=NDATA 
DO  20  I=1,N0ATA 
20  TQTAL=TOTAL+CATa( I ) 
TOTAL  =  TOFAL/RGATA 
DO  21  I=l,NDArA 
DATA( I )=OATA(I )-TOTAL 
CO  4  K=NST,NFN 
VALY=0. 
VALX=0. 
NAN=M+K+1 
NNN=K-M-1 

T AX=CAS TX»{ DATA ( NAN )+DATA{NNN)) 
TAY=CASTY»( DATA { NAN )+DATA(NNN) ) 
DO  5  J=1,M 
NAN=K+J 
NNN=K-J 

VALY=VALY+8Y( J)»(DATA(NAN)+LATA{NNN) ) 
VALX  =  VALX-i-BX(  J)*(DATA{NAN)+OATA(NNN)  ) 
VALY=(Ta«VALY+TAY)*AY 
VALX=T8«VALX+T4»DATA(K) 
VALX={ VALX+TAX)*AX 
AXT=(VALY/VALX) 
OTAN=ATANF(AXT) 

RAN(K)=SQRTF{  VALX»VALX+ VAL Y* VALY ) 
TAN(K)=SIGNF(OTAN,AXT) 


100 


21 
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C 
C 

c 
c 
c 


WRITE  OUTPUT  TAPE  6 , 6, K , DATA { K) , VALX , VALY , RAN ( K ) , TAN ( K 

) 
T(X,I5, (5(2X,E20.8)  )  ) 

7 


FORNA 
GO  TO  7 

RAN=PHASER  LENGTH 
TAN=PHASES  PROOUCtD  BY 
DUM=DUMMY  VARIAbLc  FOR 


THE  FILTER 
PLOT  DIMtNSION 


FORMAK  IHP,  14) 

CALL  POUT  (TAN,NDATA) 

CALL  DIFTN( TAN.NDATA, DA F A , RAN.BX, DUM ,BY ) 

CALL  EXIT 

END 


100 


c 
c 


ORMSBY  BAND  PASS  IN  PHASE  AND  IN  QUADRATURE  FILTERS 


C 

c 
c 
c 
c 
c 


c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


DIMENSION 
DIMENSION 
DIMENSION 
RRD^IHO 
P=3. 14159 


A( 10),H( 100) ,Z(4CCC),DUM(1C0) , OUT {4000) 
HG( ICG) ,B0UT(4C00) 
FMT  (12),  TITLE  (7) 


FMT=FORMAT  CF  DATA 
TITLE=TITLE  OF  DATA 
NDATA=NUMCER  OF  DATA 
Z=DATA  POINT 


5,25, (FMT(  I  ) ,1  =  1,12) 
5,26,(T1TLE(I ) ,1=1,7) 
5, 100,NDATA 
5,FMT, (Z(I) ,I=I,NDATA) 


READ  INPUT  TAPE 

25  FOR^'AT  (12A6) 
READ  INPUT  TAPE 

26  FORMAT  (7A6) 
READ  INPUT  TAPE 

100  F0RMAT(I4) 

READ  INPUT  TAPE 

TOTAL=0. 

RDATA=NDATA 

CC  20  I=1,NDATA 

20  T0TAL=rCTAL+Z( I) 
TOTAL=TCTAL/RDATA 
DO  21  I=1,NDATA 

21  Z(  I  )  =  Z(1)-T0TAL 


RID=LETTER  0  IN  COLUMN  I  TO  IDENTIFY  ORKSBY'S 

PARAMETER  CARD 

XN1=NUMBER  CF  TLRMS  IN  TRIGONOMETRIC  EXPANSION  OF 

TRANSFER  FUNCTION 

TD=MAXIMUM  WAVELENGTH  IN  DATA  POINTS  TO  PASS  THE 

FILTER 

TC=UPPER  CUTTOFF  WAVELENGTH  IN  DATA  POINTS 

TO=C£NT£R  WAVELLNGTH  IN  DATA  POINTS 

RK=DECIMAL  EXTENT  OF  THE  PARABOLIC  TRANSITION 

NST=CEFINES  STARTING  DATA  POINT  NUMBER  FOR  ABOVE 

PARAMETERS 

NFN=DEFINES  END  DATA  POINT  NUMBER  FOR  ABOVE  PARAMETERS 

7  READ  INPUT  TAPE  5 ,  1  ,RID , XNl , TD , TC , TO ,RK ,NST,NFN 
1  F0RMAT(A1,F9.C,4F10. 0,2110) 
IF(RID-RRD)  43,44,43 

43  WRITE  OUTPUT  TAPE  6,45 

45  FORMAT! 10X,15H  ERROR  IN  NDATA) 
CALL  EXIT 

44  IF(NST+NFN)  41,^2,41 
42  NST=1 

NFN=NDATA 
41  IF(XNl)  49,49,40 

40  WRITE  OUTPUT  TAPE  6 , 27 , ( T I TLE ( I )  ,  1  =  1  ,  7 ) 
27  FORMAT  (  1H17A6 , 6X1 3H0RMSBY  FILTER//) 

WRITE  OUTPUT  TAPE  6  ,  28 , XNl , TD, TC , TO , RK , NST ,NFN 
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28  FORMAT 
$6H  TO  = 
$  6H  RK 

WRITE  0 

29  FORMAT 
$7HVALU£ 
$  9hRESU 

NTO=TO/ 

10  T0=1.0/ 

11  TC=1.0/ 

12  TD=1.0/ 

13  TC=TO-T 

14  TD=TO-T 

15  TD=TC-T 

16  P2=P+P 

17  N1=XN1+ 

18  N2=N1+N 

19  TR=rc 

DO  2  K  = 
N=K-N1 
IF(N)  3 

3  R=N«NTO 
PN2=P2» 
RMl=1.0 
A1=SINF 
A2=SINF 
A3=SINF 
A4  =  SI;-4F 
A5=RK»R 
H  {  K  )  =  (  A 
P2NTG=P 
H{K)=H( 
HQ(K)=H 
H(K)=H( 
GO  TO  2 

4  H(K)=TC 
HQ(K)=0 

2  CONTINU 
DO  6  K  = 
T=K 
RI=0. 
RQ=C.O 
DO  5  1  = 
L=K+(I- 
IF(L)  5 

8  IFtNCAT 

9  RI=RI+H 
RG=RG+H 

5  CONTINU 
OUT(K)= 
BOUT(K) 
WRI  TE  0 

22  FORMAT! 

6  CONTINU 


(7H  XNl  =  F10.6,6H  TO  =  F10.6,6H  TC  =  F10.6, 

F10.6, 
=  F1C.6  /  7H  iMST  =  16, 7H  KFN  =  lb///) 
UTPUr  TAPE  6,29 
{X6HNUMBhR6X10HlNPUT  CATA13X7hV ALUE  Xl^X 

Y14X 
LTANTlSXbHTHETA//) 
(XNl+XNl) 
TO 
TC 
TD 
C 
D 
C 


1.0 
1-1 

1,N2 

»4,3 

R 

-RK 

(PN2*(TC+RM1*TR) ) 

( PN2*(TC+RK»TR) ) 

(PN2»{TC+TK) ) 

(PN2*TC) 

M1»TR«TR»P»R»PN2*PN2 

1+A2-A3-A4)/A5 

N2»TC 

K)+h(K) 

(K)»SINF(P2NT0) 

K)«CCSF(P2.NTU) 

+TC+TR 


iMST.NFN 


1,N2 

N1)«N 

,5,8 

A-L) 

(I)»Z 

Q(  I)» 

c 

SQRTF 

=  ATAN 

UTPUT 

X,I5, 

E 


TO 

5,5,9 

(L) 

Z(L) 

(RI*RI+RQ»RQ) 
F(RQ/RI ) 

TAPE  6,22,K,Z(K)  ,RI  ,RQ,CUT(K) ,BCUT(K) 
(5(2X,E20.8)  )  ) 
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C 

c 
c 
c 
c 


GO  TO  7 

UUT=PHASER  LENGTH 

BOUT=PHASES  PKODUCEO  BY  FILTER 

DUM=CUMMY  VARIABLE  FOR  PLOT  DIMENSION  IOC 

49  CALL  RPLTIOOT.NDATA.DUM) 
CALL  RPLTIBCUT.NDATA.DUM) 
WRITE  OUTPUT  TAPE  6  ,  30  ,  ( TI TLE { I )  ,  I  =  1  ,  7  ) 

30  FORMAT  (1HP7A6 , 6X13H0RMSBY  FILTER) 
WRITE  OUTPUT  TAPE  6,31,NCATA 

31  FORMAT  (1HPI4) 

CALL  POUT(BCUT,NDATA) 

CALL  DIFTN(B0UT,NDATA,2,0UT,H,DUM,HQ) 

CALL  EXIT 

END 
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C  PROGRAM  DIFTN  USING  LANCZOS  FGRNULA  FOR 

C  DIFFERENTIATION  KITH  WEIGHTS  CHANGED  TO  ZERU   BEYOND 

C  END  POINTS 

C  PROGRAM  CORRECTS  FOR  CYCLIC  DISCONTINUITY  OF  THE  ARC 

C  TANGENT  FUNCTION 

C 

1  F0RMAT(I4) 

2  F0RMAT(6£12.0) 

3  F0R^'AT(5E20.5) 
6  F0RMAT{2F10.C) 

Al  FORMAT! IHl,//, IGX,26H  RDATA , FRING , W ICTH ,DX ,CCF  =  ,  5E  1  8  .5 

) 

42  F0RFAT(12A6) 

43  F0RMAT(/,10X,12A6) 

40  FORMAT!///, I8X,2H  X,I6X,4H  TAN,15X,5h  TAN2,15X,5H  DTAN 

t 
$I2X,8H  STRAINS) 

DIMENSION  DrANl4C00) ,DTR(1CG) 

DIMENSION  TAN(4G00),TAN2{4CGC),A{100),TR(100) 
C 

C      KR=NUMBER  OF  POINTS  IN  POLYNOMIAL  APPROXIMATION  IN  THE 
C      NEIGHBORHOOD  OF  THE  POINT  WHERE  THE  DERIVATIVE  IS 
C      SOUGHT 

C      FRING=NUMBER  OF  LINES  PER  INCH  OF  THE  MASTER 
C      WIDTH=LENGTH  OF  THE  MODEL  IN  THE  DIRECTION  OF  SCAN 
C      A=TITLE  OF  THE  DATA 

C      NDATA=RUATA=NUMBER  OF  DATA  POINTS 
C      TAN=PHASES  GENERATED  BY  FILTER 
C 

READ  INPUT  TAPE  5,I,KR 

RcAC  INPUT  TAPE  5 , 6, FR ING, W I D TH 

READ  INPUT  TAPE  5, 42 , ( A ( I ) ,  1  =  1 , 1 2 ) 

READ  INPUT  TAPE  5,1,NDATA 

READ  INPUT  TAPE  5 , 2, ( TAN ( I ) , I =1 , NDATA ) 

PI=3.1416 

NT=NDATA/100 

RDATA^NDATA 

DX=WIDTH/RDATA 

CCF=I.0/{2.C»FRING»PI*DX) 

WRITE  OUTPUT  TAPE  6 , 41 , RCATA , FRING, W IDTH , DX , COF 

WRITE  OUTPUT  TAPE  6 , 43  ,  { A (  I  )  ,  1  =  I , 12 ) 

WRITE  OUTPUT  TAPE  6,40 

AMUL=0. 

DO  5  K=l, NDATA 

A1=TAN{K)-TAN(K-1) 

AI=SIGNF(2.,A1) 

A2=TAN{K-l)-TAN(K-2) 

A2=SIGNF( l.,A2) 

A3=SIGNF(4. ,TAN(K) ) 

A4  =  SIGNF(3.  ,TAN(K-I)  ) 

AA=(A1+A2+A3+A4) 

AA2=ABSF{AA) 

IF(AA2-3.0)  10,11,11 
10  IF(AA)  12,11,13 


lOli 


c 
c 
c 
c 
c 

c 
c 
c 
c 
c 
c 


12  AMUL=AMUL+PI 
GO    TC    11 

13  AMUL=AMUL-PI 

11  TAN2(K)=TAN(K)+AMUL 
5  CONTINUt 
TERr'^O. 
DO  20  1=1, KR 
ALPHA=I 

20  TERM=TERM+ALPHA«ALPHA 
TERM=TERM»2.0 

K  =  0 

X  =  0. 

DO  21  I=1,NDATA 

R=0. 

DO  22  J=1,KR 

ALPhA=J 

IPJ=I+J 

If<J=I-J 

IF(IMJ)  22,33,33 
33  ITEM=NCATA-IPJ 

IF(  ITEM)  22,32,32 
32  R=R+ALPHA»TAN2( i P J )-ALPHA»TAN2 ( I FJ ) 
22  CONTINUE 

X=X+CX 

STAN=R/TERM 

DTANl I )=STAN*COF 

TAN2=PHAStS  CORRECTED  FOR  Dl SCONTINUl TI tS 

STAN=UNCORRECTEtj  STRAINS 

DTAN=STRAINS 

21  WRITE  OUTPUT  TAPE  6 , 3, X , TAN { I ) , TAN2 ( I ) , STAN , DTAN ( I) 

DTR=CTAN  REDUCED  IN  NUMBER  FOR  PLOTTING,  ALWAYS  100 

VALUES 

TR=TAN2  REDUCED  AS  ABOVE 

A=DUNMY  VARIABLE  FOR  PLOTTING  DIMENSION  100 

J  =  0 

CO  A  1=1, IOC 
J=J+NT 

DTR( I  )  =  DTAN( J) 
4  TR{  I  )  =  TAN2( J) 

CALL  RPLT(TR,  100,A) 
CALL  RPLT(DTR,1G0,A) 
RETURN 
END 
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c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


c 
c 
c 
c 
c 
c 
c 


SUBROUTINE  DIFTN  USING  LANCZOS  FCRKULA  FOR 

DIFFERENTIATION  WITH  FUNCTION  REFLECTED  K.R.T.  END 

POINTS 

PROGRAM  CORRECTS  FOR  CYCLIC  DISCONTINUITY  OF  THE  ARC 

TANGENT  FUNCTION 

TAN=PHAS£S  GENERATED  BY  FILTER 

NDATA=RDATA=NUMBER  OF  DATA  POINTS 

DTAN=STRAINS 

TAN2=PHASES  CORRECTED  FOR  DISCONTINUITIES 

CTR=DTAN  REDUCED  IN  NUMBER  FOR  PLOTTING,  AL^^AYS  100 

VALUES 

TR=TAN2  REDUCED  AS  ABOVE 

A=DUMMY  VARIABLE  FOR  PLOTTING  DIMENSION  100 

SUBROUTINE  D  IFTN ( TAN , NDATA , DTAN , TAN2 ,DTR , A , TR) 

1  FORPAT(M) 

2  FCRMAT(6E12.Q) 

3  F0RMAT(!3E20.5) 
6  FORMAT(2FI0.0) 

41  FORNATC 1H1,//,I0X,26H  RDATA  ,  FRING , WI CTH , DX ,CGF  =  , 5E 18 . 5 

) 

42  FORMAT! 12A6) 

43  FORMAT(/,10X,12A6) 

40  FOR^ATi///, 18X,2H  X,16X,4H  TAN,15X,5H  TAN2,15X,5H  DTAN 

t 
$12X,8H  STRAINS) 
DIMENSION  DTAN{4000),DTR( 100) 
DIMENSION  TAN(4000),TAN2(4COC),A(100),TR(100) 

KR=NUMBER  OF  POINTS  IN  POLYNOMIAL  APPROXIMATION  IN  THE 

NEIGHBORHOOD  OF  THE  POINT  WHERE  THE  DERIVATIVE  IS 

SOUGHT 

FRIISG=NUMBER  OF  LINES  PER  INCH  OF  THE  MASTER 

WIDTH=LENGTH  OF  THE  MODEL  IN  THE  DIRECTION  OF  SCAN 


READ 
READ 
PI  =  3 
NT=N 
RDAT 
DX  =  W 
COF  = 
WRIT 
WRIT 
AMUL 
DO  5 
A1  =  T 
A1  =  S 
A2  =  T 
A2=S 
A3  =  S 
A4  =  S 
AA=( 


INPU 
INPU 
.1416 
DATA/ 
A=NCA 
IDTH/ 
1.0/( 
E  CUT 
E  OUT 
=  0. 

K=l, 
,AN(K) 
IGNF( 
AN(K- 
IGNF( 
IGNF( 
IGNF{ 
A1  +  A2 


T  TAPE  5,1,KR 

T  TAPE  5, 6, FRING, WIDTH 


100 

TA 

RDATA 

2.C*FRING*PI*DX) 

6,41,RDATA,FRING,WIDTH,DX,C0F 

-  -in 


PUT 
PUT 


TAPE 
TAPE 


6,40 


NDATA 
-TAN(K-1} 
2., AD 
l)-TAN(K-2) 
l.,A2) 
4.,TAN(K)  ) 
3.,TAN(K-1)) 
+A3+A4) 
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c 
c 
c 


AA2  =  /kBSF(AA) 
IF(AA2-3.0)     10,11,11 

10  IF(AA)     12,11,13 

12  AMUL=AMUL+PI 
GO    TO     11 

13  AMUL=AMUL-PI 

11  TAN2(K)=TAN(K)+AMUL 
5  CONTINUE 

TERM=0. 

DO  20  1=1, KR 

ALPHA=I 

20  TERM=TERM+ALPHA*ALPHA 
TERM=TERM*2.0 

X  =  0. 

DO  21  I=1,NDATA 

R  =  0. 

DO  22  J=1,KR 

ALPHA=J 

IPJ=I+J 

IMJ=I-J 

IFl  IMJ)  31, 33,33 

31  IMJ^-IKJ 

33  ITEM=NDATA-IPJ 

IF(ITEM)  32,22,22 

32  IPJ=NDATA+ITEM 

22  R=R+ALPHA*TAN2( IP J )-ALPHA»TAN2 { IMJ) 
X=X+CX 
STAN=R/T£RM 
DTAN{ I )=STAN»CGF 

STAN=UNCaRRECTEL-  STRAINS 

21  WRITE  OUTPUT  TAPE  6, 3, X , TAN ( I ) , TAN2 ( I ) , STAN , DTAN ( I ) 
J  =  0 

DO  A  1=1, ICC 
J=J*NT 

OTR( I )=DTAN( J) 
4  TR(I)=TAN2(J) 

CALL  RPLT(TR,10U,A) 

CALL  RPLT(0TR,1U0,A) 

CALL  EXIT 

RETURN 

END 
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C      SUBROUTINE  DIFTN  USING  LANCZQS  FORMULA  FCR 

C      DIFFERENTIATION  KITH  WEIGHTS  CHANGED  TO  ZERO  BEYOND 

C      END  POINTS 

C      PROGRAM  CORRECTS  FOR  CYCLIC  DISCONTINUITY  OF  THE  ARC 

C      TANGENT  FUNCTION 

C 

C      TAN=PHASES  GENERATED  BY  FILTER 

C      NDATA=RDATA=NUMBER  OF  DATA  POINTS 

C      DTAN=STRA1NS 

C      TAN2=PHASES  CORRECTED  FOR  DISCONTINUITIES 

C      CTR=DTAN  RECUCEC  IN  NUMBER  FUR  PLOTTING,  ALWAYS  100 

C      VALUES 

C      TR=TAN2  REDUCED  AS  ABOVE 

C      A  =  DUMMY  VARIABLE  FOR  PLUTTI.NG  DIMENSION  100 

C 

SUBROUTINE  DIFTN(TAN,NDATA,DTAN,TAN2,0TR,A,TR) 

1  FORMAT (14) 

2  FCRMAT(6E12.0) 

3  FORMAT(5E20.5) 
6  FORMAT(2F10.0) 

41  FORMAT{ 1H1,//,1GX,26H  RDATA , FRING, W I CTH ,DX ,CGF=, 5E 18. 5 

) 

42  FORMAT( 12A6) 

43  FORMAT!/, lOX, 12A6) 

40  FORMAT(///, 18X,2H  X,16X,4H  TAN,15X,5H  TAN2,15X,5H  DTAN 

t 
$12X,BH  STRAINS) 

DIMENSION  DTAN(4000),DTR(100) 

DIMENSION  TAN(4C00),TAN2(4000),A(1C0),TR(100) 
C 

C      KR=NUMBER  OF  POINTS  IN  POLYNOMIAL  APPROXIMATION  IN  THE 
C      NEIGHBORHOOD  OF  THE  POINT  WHERE  THE  DERIVATIVE  IS 
C      SOUGHT 

C      FRING=NUMBER  OF  LINES  PER  INCH  OF  THE  MASTER 
C      WIDTH=LENGTh  OF  THE  MODEL  IN  THE  DIRECTION  UF  SCAN 
C 

REAC  INPUT  lAPE  5,1,KR 

REAC  INPUT  TAPE  5 , 6 , FR ING, W IDTH 

PI=3.14I6 

NT=NCATA/100 

RDATA=NCATA 

DX=WIDTH/RDATA 

C0F=1.0/{2.C»FRING»PI*DX) 

WRITE  OUTPUT  TAPE  6 , 41 , RDATA , FRING, W IDTH , OX ,COF 

WRITE  OUTPUT  TAPE  6,40 

AMUL=0. 

CO  5  K=1,NDATA 

A1=TAN(K)-TAN(K-1) 

A1=SIGNF(2. ,A1) 

A2=TAN(K-l)-TAN(K-2) 

A2=SIGNF( l.,A2) 

A3  =  SIGNF14.,TAN(K)  ) 

A4=SIGNF{3. ,TAN(K-1)  ) 

AA=(A1+A2+A3+A4) 
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c 
c 
c 


AA2=A8SF(AA) 
IF(AA2-3.0)     10,11,11 

10  IF(AA)     12,11,13 

12  AMUL=AMUL+PI 
GO    TC    11 

13  AKUL^AMUL-PI 

11  TAI\i2(K)  =  TAN(K)+AMUL 
5  CONTINUE 

TERM=0. 

DO  20  1=1, KR 

ALPHA=I 

20  TERM=TERM+ALPHA*ALPHA 
TERM=rERM»2.0 

X  =  0. 

CO  21  I  =  l,fMDATA 

R  =  0. 

CO  22    J=1,KR 

ALPHA=J 

IPJ=I+J 

IMJ=I-J 

IF(If'J)  22,33,33 
33  ITEf^  =  NCATA-IPJ 

IF(ITEM)  22,32,32 
32  R  =  R  +  ALPhA»TAN2( 1 P J ) -ALPHA»TAN2  (I  M J  ) 
22  CONTINUE 

X=X+DX 

STAN=R/TERM 

CTAN( I )=STAN»COF 

STAN=UNCORRECTED  STRAINS 

21  WRITE  OUTPUT  TAPE  6, 3, X , TAN ( I ) , TAN2 ( I ) , STAN ,OTAN ( I ) 
J  =  0 

DO  4  1=1,100 
J=J+NT 

CTR( I)=OTAN(J) 
4  TR{  I)  =  TAN2( J) 

CALL  RPLT(TR,10G,A) 
CALL  RPLT(DTR, 100, A) 
RETURN 
END 
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C      SUBROUTINE  PGUT  PUNCHES  DATA  GENERATED  BY  PROGRAMS 

C-     POINT=CATA  TC  Bt  PLOTTED 

C      N  =  NUf^BER  OF  DATA  POINTS  TO  BE  PLOTTED 

C 

SUBROUTINE  POUT ( PO INT, N ) 

DIMENSION  POINTCI) 

J  =  N/6 

JG=1 

DO  I  1=1, J 

JE=J0+5 

WRITE  OUTPUT  TAPE  6 , 2  ,  ( POI NT ( JT ) , JT  =  JO , JE ) 

1  J0=JE+1 

WRITE  OUTPUT  TAPE  6 , 2 , ( PCI  NT ( JT ) , JT= JO , N ) 

2  FORMAT( 1HP,6E12.5) 
RETURN 

END 


no 


c 
c 
c 
c 
c 


SUBRCUTiNt  RPLT  PLOTS  DATA  GENERATED  BY  PROGRAMS 

POINT=DATA  TO  Bt  PLOTTED 

N=NUMBER  OF  DATA  POINTS  TO  BE  PLOTTED 

A  =  DU^'MY  VARIABLE  FOR  PLOT  DIMENSION  100 

SUBROUTINE  RPLT  {  PO  IiN  T  ,  N  ,  A  ) 
DIMENSION  POINT(i) ,A(1) 
B=1H 
X=1HX 

CO  6  1=1,100 
6  A( I )=B 

GP=PCINT( 1) 
GN=PCINT(1) 
CO  2  1  =  1, N 
IF(POINT(I)-GP)  3,2,4 

4  GP=POINT(I) 
GO  TC  2 

3  IF(PCINT( I )-GN)  5,2,2 

5  GN=PCINT(I) 
2  CONTINUE 

R=998./(GP-GN) 

NS=-GN*R+1. 

WRITE  OUTPUT  TAPE  6, 8 , R , GP , GN ,NS 

8  FORMAT(  lOHlRATIb  I  S  ,  E  14  .  7,  1  31-   MAXIMUM  IS,E14.7, 
$13H   MINIMUM  IS 

lfcl4.7,17h   SCALED  ZERO  IS, 14,//) 
DO  9  1  =  1, N 

IOT=(PCINT( I )-GN)*R+1.0 
J=( ICT/1C)+1 
A(J)=X 

WRITE  OUTPUT  TAPE  6, 1 1 , I CT , ( A ( K )  , K=i , I  00 ) 
11  FORMATt IX, I4,5X, ICOAl) 

9  A(J)=B 
RETURN 
END 
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Fig.   1    Modulation  Transfer  Function  Measurement 
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SENSITOMETRIC 
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OF  A  NEGATIVE  ASSUMING 
D  =  LOG  E 


TRANSMITTANCE  VARIATION  OF  THE 
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Fig.  3  Spatial  Waveform  Distortion  Due  to  the  Film 
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Fig.  10  Coherent  and  Incoherent  Light  Systems 
a)  Coherent,  b)  Coordinates,  c)  Incoherent 


122 


Fig.  11  Coherent  Optical  Filtering  by  Removal  of  the  Bias  Term 
of  the  Grid.   (Insert  above:   50X  magnification  of  the 
master  grid  lines  before  filtering.  Insert  below: 
$0X  magnification  of  the  master  grid  lines  after 
filtering.) 
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Fig.  17  The  Electro -Optical  System 
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Fig.  18  Ring  Under  Diametral  Compression 
(Moire  of  U) 
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Fig.  hi     Total  Residual  Strains  in  Ring  No,   2  by  the  Moire  Method 
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Fig,  Ii2     Boring-Out  Test  Set-Up 
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Figo  h5     Residual  Plastic  Strains  in  Ring  No.  2 
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Fig,  16   lOOX  Magnification  of  an  Attempt  to  Engrave  750  Lines 
per  Inch  on  an  Insufficiently  Polished  Aluminum  Sample 
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Fig,  li7   hOX  Magnification  of  a  Photoengraved  Cross  Grid  of  300 
Lines  per  Inch 
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Fig,   Ii9       )iOX  Magnification  of  an  Etched  Cross  Grid  of  300  Lines 
per  Inch  After  Removal  of  the  Resist 


APPENDIX  E 
DETERMINATION  OF  FILTER  PARAMETERS 

Numerous  comment  cards  have  been  included  in  the  programs  of 
Appendix  C  to  aid  the  user  and  to  clarify  the  meaning  of  the  abbrevia- 
tions of  the  Fortran  Statements. 

We  will  describe  here  the  sequence  of  operations  leading  to 
the  determination  of  the  adequate  values  for  the  independent  parameters 
of  Filters  1  and  2. 

Filter  1 

This  filter  possesses  two  independent  parameters,  the  center 
frequency  f  and  the  number  of  terms  in  the  trigonometric  approximation 
of  the  transfer  function  m. 

It  is  important  to  notice.  Figure  21,  that  the  band  width, 
Af  =  2TT/m,  and  the  roll-off  frequency  interval,  2Tr/ra,  depend  on  the 
value  of  m. 

.  Keeping  this  in  mind  we  proceed  as  follows: 

a)  Display  in  a  standsird  IBM  80  x  80  listing  the  digitalized 
values  of  the  light  intensity  trace  obtained  at  equally  spaced  intervals 
along  the  direction  of  scan. 

b)  DeteiTnine  the  shortest  and  longest  wavelength  of  the  fringes 
in  the  trace  in  data  point  numbers.  The  difference  of  these  two  quanti- 
ties defines  the  wavelength  deviation  AA.  The  frequency  deviation  is 
obtained  as  Af  =  2n/AX. 
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c)  Obtain  the  center  wavelength  X  of  the  trace  as  the  average 
of  the  shortest  and  longest  wavelengths  determined  in  (b) .  The  center 
frequency  will  then  be  f  =  2n/X 

d)  Assume  for  the  first  trial  run  that  the  band  width  of  the 
filter  is  twice  the  frequency  deviation.  Then 

A.    m 


and  the  parameter  m  can  be  defined  by 

TT    A\ 

"•  =  Af  =  T 


(E.2) 


e)  Load  the  calculated  parameters  f  and  m  into  the  program  of 
Page  97  and  perform  the  first  trial  run. 

f )  Inspect  the  plotted  output  of  the  amplitude  vector  that  has 
been  called  phaser  length  in  the  program.  These  valued  should  be  approx- 
imately constant  for  minimum  error  in  a  properly  tuned  filter.  Inspect 
the  plotted  output  of  the  displacement  curve  to  determine  the  wave- 
length of  the  noise,  if  any,  that  is  passing  through  the  filter. 

g)  If  unsatisfactory  results  have  been  obtained  in  the  first 
trial,  adjust  the  values  of  f  and  m  in  successive  runs  until  the 
filter  yields  a  smooth  displacement  curve  with  minimum  variations  in 
the  magnitude  of  the  amplitude  vector. 

In  evaluating  this  filter  for  a  possible  moire  application  the 
reader  should  keep  in  mind  the  following: 

Some  sfjurious  disturbances  in  the  output  of  the  filter  can  be 
due  to  noise  content  close  to  the  signal  frequencies  that  are  passing 
through  the  filter  in  the  roll-off  interval  2TT/m.  In  order  to  narrow 
this  interval  m  can  be  Increased.  However,  this  decreases  the  band 
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width  and  might  obliterate  some  information  on  the  desired  signal. 

On  the  other  hand,  if  the  useful  signal  is  contained  in  a  broad 
band,  the  value  of  ra  necessarily  must  be  small  to  widen  the  band  pass 
of  the  filter.  This  will  cause  the  roll-off  interval  to  widen  and  will 
reduce  the  number  of  terms  in  the  series  approximation  of  the  transfer 
function  imparing  the  accuracy  of  the  procedure. 

The  shortcomings  mentioned  above  lead  us  to  believe  that  the  use 
of  Filter  1  should  be  restricted  to  moire  light  intensity  traces  of 
narrow-band  frequency  composition  with  noise  spectra  well  removed  from 
the  signal  spectra. 

Filter  2 

This  filter  possesses  five  independent  parameters,  the  number 
of  terms  in  the  trigonometric  approximation  of  the  transfer  function  N, 
the  center  frequency  f  ,  the  cut-off  frequency  f  ,  the  termination  fre- 
quency f  and  the  coefficient  k  that  defines  the  extent  of  the  parabolic 
transitions  of  the  roll-off,  Figure  22. 

This  filter  affords  the  possibility  of  defining  the  number  of 
terms  in  the  expansion  of  the  transfer  function,  the  band  width  and  the 
roll-off  frequency  interval  independently. 

The  independent  parameters  are  determined  as  follows: 

a),  b)  and  c)  are  the  same  as  for  Filter  1. 

d)  Assume  for  the  first  trial  run  that  the  cut-off  frequency  is 


^0  =  ^o  *  T  =  ^o  ^  s;  (E-3) 


e)  The  values  for  N  and  k  will  depend  on  the  noise  composition 
of  the  input;  however,  N  will  rarely  exceed  I;0  and  a  value  of  k  between 
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0.2  and  O.U  will  produce  good  results  in  most  cases. 

f)  The  parameter 

f  -  f 

s 

that  permits  us  to  define  the  termination  frequency  f ,  as  a  function 

of  the  cut-off  frequency  f  and  the  sampling  frequency  f  ,  can  be 

c  s 

obtained  by  using  the  criteria  given  in  [30]. 

g)  Load  the  calculated  parameters  f  ,  N,  k,  f  and  f  into  the 
program  of  Page  100  and  perform  the  first  trial  run. 

h)  Perform  the  inspections  described  in  (d)  of  Filter  1. 

i)  If  unsatisfactory  results  are  obtained  in  the  first  trial, 
adjust  the  values  of  the  independent  parameters  in  successive  runs 
until  the  filter  is  properly  tuned. 

The  existence  of  five  independent  parameters  will  permit  more 
accurate  and  fast  tuning  of  this  filter  than  of  Filter  1. 
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